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1. Executive summary 



The phenomenon of neutrino oscihations, arguably the most significant advance in particle physics 
over the past decade, has been established through measurements on neutrinos and anti-neutrinos 
produced in the sun, by cosmic-ray interactions in the atmosphere, nuclear reactors, and beams 
produced by high-energy particle accelerators. In consequence, we know that the Standard Model is 
incomplete and must be extended to include neutrino mass, mixing among the three neutrino flavours, 
and therefore lepton-flavour non conservation. These observations have profound implications for the 
ultimate theory of particle interactions and for the description of the structure and evolution of the 
Universe. 

These exciting possibilities justify an energetic and far reaching programme, an essential part of 
which is to make precision measurements of the oscillation parameters. Assuming the three flavours 
and the unitary neutrino-mixing matrix that is presently favoured, oscillation measurements can be 
used to determine the three mixing angles and the critical phase parameter that can provide a new 
source of CP-invariance violation. Neutrino oscillation measurements can also be used to determine 
the two (signed) mass differences. This programme is similar to the long-standing investigations of 
quark mixing via the CKM matrix and it would now seem to be clear that an understanding of the 
flavour problem will definitely need precision measurements in both quark and lepton sectors. 

Not all the properties of the neutrino can be determined by oscillation experiments. Equally impor- 
tant is the determination of the Major ana or Dirac nature of the neutrino which requires the ongoing 
and planned neutrinoless double beta decay experiments. In addition, although oscillation measure- 
ments determine the mass differences, they are insensitive to the absolute mass, mi, of the lightest 
mass state. The determination of mi requires a very precise measurement of the end-point of the 
electron spectrum in beta decay. 



Coordination and Outreach 

The members of WPG have held meetings during the general meeting at Strasbourg in June 2010, 
at the CERN workshop NuThemes in September 2010 and at the Rutherford lab during the IDS- 
NF meeting also in September 2010. A summary of the EUR0NU-WP6 workshop held at from 810 
June 2009 at Cosener's House, Abingdon, UK on the subject "Flavour physics in the era of precision 
neutrino experiments" has been included in this report. 

The group has kept close contact with the IDS-NF concerning the physics of the neutrino factory, 
with WP4 as regards the scenarios for the beta-beam and with WPS as regards detector performance 
and systematic errors. During 2010, the WP6 results of the first year of the project (2009) were 
summarized in a report that was submitted to the archives to inform the wider community [20| . 
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Global analysis and interpretation of present data 

The members of WPG have carried out a number of important studies for the EURONU project. 
The first is a revision of the global analysis of neutrino oscillation data including the latest results, 
such as those from MiniBOONE and MINOS, as well as new cosmological data. The values of the 
atmospheric and solar parameters and the bounds on the angle ^13 are of course essential for the 
physics optimization of future facilities. Obviously if there is physics beyond the standard three- 
neutrino scenario, for example the presence of new sterile species, this could also change drastically 
the optimization of baseline and energy of the neutrino beam. 

Evaluation of physics performance, optimization and compcirison 

A number of studies of the physics potential of the future facilities as regards the standard three- 
neutrino oscillation scenario, as well as new physics, have been performed during the past year. The 
results have been published in peer-review journals and are also listed as EURONU documents, or 
will be presented here as internal documents. Between these, the report contains a new design of the 
SPL-Frejus Super-Beam is included (see Ref. 0] for more details), the study of the potential of high-7 
electron-capture beta-beams in the context of the LAGUNA european project has also been reviewed 
0], a review of all beta-beam setups that have been discussed in the literature in the past 0] and a 
dedicated study of atmospheric neutrino backgrounds at the ICAL@INO detector exposed to a high-7 
/3-beam 

Tools for physics studies 

A new release of the GLoBES package including migration matrices for signal and background was 
made available during 2010. 
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2. The Physics case 



The main motivation of a future neutrino physics programme is to unveil what the new physics 
associated to neutrino masses is. We know for sure that new degrees of freedom must be added to 
the Standard Model (e.g. right-handed neutrinos) at some energy scale A. If A is much larger than 
the electroweak scale, there is a natural explanation of why neutrinos are so light. Indeed the effects 
of any such new physics must be generically well described at low energies by an effective Lagrangian 
which contains the Standard Model, plus a tower of higher dimensional operators constructed with 
the SM fields and satisfying all the gauge symmetries: 



The effective operators, Oj, are ordered by their mass dimension, since the higher the dimension, the 
higher the power of A that suppresses them. The dominant operator is therefore the lowest dimensional 
one, with d = 5, which is precisely the Weinberg's operator: 



which, as is well known, induces three new ingredients to the minimal SM: 

• Neutrino masses 

• Lepton mixing 

• Lepton number violation 

In this context, neutrino masses are very small, because they come from an effective operator which is 
suppressed by a high energy scale. If we go to operators of d = 6, that are suppressed by two powers 
of A, these will generically induce new physics in dipole moments, rare decays, etc. Beyond d = 6 we 
would find operators inducing non-standard neutrino interactions (NSI). 

It is also possible that the scale A is at or below the electroweak scale, or in other words that neutrino 
masses are linked to light degrees of freedom, i.e. a hidden sector which we have not detected yet, 
because it is weakly interacting. Such scenarios do not offer an explanation of why neutrinos are light, 
but neutrinos are the natural messengers with such hidden sectors, since they are the only particles 
in the SM carrying no conserved charge. Such new physics could be related to other fundamental 
problems in particle physics such as the origin of dark matter and dark energy. 

Even though it is not guaranteed that we can fully understand the new physics associated to neutrino 
masses by measuring them, it is quite clear that we have a good chance to learn something more about 
it by testing the Standard scenario of mixing with future and more precise neutrino experiments. 
In particular we should be able to measure all the fundamental parameters: three mass eigenstates 
(mf,m|,m|), three angles (^12, ^13) ^23) and one or three CP-violating phases {6,ai,a2)- But, also, it 
will be very important to search for new physics beyond neutrino masses and mixings, in particular 
for those effects that are generic in many models of neutrino masses, such as violations of unitarity, 
non-standard interactions or the presence of light sterile species. To some extent these searches can 
also be improved in future facilities and this should be evaluated. Typically such analyses imply 
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dealing with a much larger parameter space, which calls for new tools to perform the fits, in particular 
Montecarlo methods. 

Many studies in the last ten years have shown that we can measure ^13 , discover leptonic CP violation 
and determine the neutrino hierarchy in more precise neutrino oscillation experiments, searching for 
the subleading channel Vf. -H- or its CP-conjugate channel ^ iii the atmospheric range. In this 
first section, we present the results obtained within the work of EUR0nu-WP6 in 2010 concerning: the 
status of leptonic mixing global fits [it; the prospects for ^13 searches [3]; and, the direct determination 
of solar neutrino fluxes from solar neutrino data 



2.1. The leptonic flavour sector pj 



It is now an established fact that neutrinos are massive and leptonic flavors are not symmetries of 
Nature 



21 



22l |. In the last decade this picture has become fully proved thanks to the upcoming of 



a set of precise experiments. In particular, the results obtained with solar |23h32|| and atmospheric 



neutrinos 



33, 
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nuclear reactors 



have been confirmed in experiments using terrestrial beams: neutrinos produced in 
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of hundreds of kilometers 



3g] and accelerators 37H40f| facilities have been detected at distances of the order 



The minimum joint description of all the neutrino data requires mixing among all the three known 



neutrinos {v^ 



which can be expressed as quantum superpositions of three massive states Ui 



{i = 1,2,3) with masses mj . This implies the presence of a leptonic mixing matrix in the weak charged 



current interactions 
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quark sector, there are two physical phases r/j associated to the Majorana character of neutrinos and 
which are not relevant for neutrino oscillations 
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Given the observed hierarchy between the solar and atmospheric mass-squared splittings there are 
two possible non-equivalent orderings for the mass eigenvalues, which are conventionally chosen as 

mi < 1712 < nT-s with Arn^i ^ l^^iii 



7713 < mi < 1712 with Am^x <ti |Am|x 



Amial 



and 
and 



Amii > ; 
Am|i < . 



(4) 
(5) 



As it is customary we refer to the first option, Eq. as the normal (N) scheme, and to the second 
one, Eq. ([5]), as the inverted (I) scheme; in this form they correspond to the two possible choices of 
the sign of Am|^. In this convention the angles 6ij can be taken without loss of generality to lie in 
the first quadrant, Oij G [0,7r/2], and the phases 6cp, iji S [0,27r]. 

Thanks to the synergy amongst a variety of experiments involving solar and atmospheric neutrinos, 
as well as man-made neutrinos at nuclear power plants and accelerators, we have now a relatively 
detailed picture of the parameters describing three-flavor neutrino oscillations 
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Figure 1. Allowed parameter regions (at 90%, 95%, 99% and 99.73% CL for 2 d.o.f.) from the combined 
analysis of solar data for ^13 ~ 0. The best-fit point is marked with a star. For comparison we also 
show as empty regions (the best-fit is marked by a circle) the results prior to the inclusion of the latest 
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32| and the low energy threshold 



Ga capture rate of SAGE [25| , the energy spectrum of Borexino 

analysis of the combined SNO phase I and phase II [sB]. In both analysis we use as inputs the GS98 



solar model fluxes and the Gallium capture cross-section of Bahcall [49 1 



2.1.1. Leading Arn^i oscillations: solar and KamLAND data 



In the analysis of solar neutrino experiments we include the total rates from the radiochemical ex- 
periments Chlorine 23(], Gallex/GNO 2J] and SAGE [25[. For real-time experiments we include the 
44 datapoints of the electron scattering (ES) Super-Kamiokande phase I (SK-I) energy-zenith spec- 



trum [26( 1 and the data from the three phases of SNO 27l429l| , including the results on the low energy 
threshold analysis of the combined SNO phase I and phase II 30(] (which we label SNO-LETA). We 
also include the main set of the 192 days of Borexino data [311] (which we label Borexino-LE) as well 
as their high-energy spectrum from 246 live days [32l | (Borexino-HE). 

In Fig. [1] we show the present determination of the leading parameters Am2i and 612 from the 
updated oscillation analysis of the solar neutrino data described above in the context of the GS98 
solar model. For comparison we also show the results obtained prior to the inclusion of the latest Ga 



capture rate of SAGE [25|], the energy spectrum of Borexino [3ll.l32l| and the SNO-LETA results [30(] for 



the same solar model. As seen in this figure, the inclusion of these results lead to an improvement on 
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Figure 2. Allowed parameter regions (at 90%, 95%, 99% and 99.73% CL for 2 d.o.f.) from the analysis 
of atmospheric data (full regions, best-fit marked with a star) and LBL data (void regions, best-fit 
marked by a circle) for 6*13 = and Am^i = 7.6 x fO~^ eV^. 



the determination of both 612 and Aml^ and for this last one the best-fit value slightly increases. The 
most quantitatively relevant new information arises from the inclusion of the SNO-LETA results. The 
inclusion of Borexino tends to shift the region towards slightly lower values of 0i2 angle. Conversely, 
if the analysis is done in the context of the AGSS09 model the region is shifted towards slightly larger 

^12- 



2.1.2. Leading Am|^ oscillations: atmospheric and accelerator data 



In this section we present two different analyses of the atmospheric data. The first one is very similar 
to the one detailed in Ref. 



4l[, and includes the results from the first run of Super-Kamiokande, w^ 



lich 



accumulated data from May 1996 to July 2001 (1489 day exposure) and is usually referred as SK-I [33l |. 
as well as the data obtained with the partial coverage after the 2001 accident (804 day exposure), the 
so-called SK-II period [50|. We will refer to this analysis as SK(I-|-II). The second one is based on 
the new analysis recently presented by the Super-Kamiokande collaboration including also the data 
taken from December 2005 to June 2007, usually referred as SK-III [3j]. Apart from the inclusion 
of these new event rates, in this data release the selection criteria and the corresponding estimate of 
uncertainties for the SK-I and SK-II periods have been changed with respect to the previous SK(I-|-II) 
analysis. We have therefore performed a reanalysis of the new combined samples from phases I, II 
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and III as presented in [3j]. We refer to the results of this analysis as SK(I+II+III). It is important 
to point out that already since SK-II the Super-Kamiokande collaboration has been presenting its 
experimental results in terms of a large number of data samples. The rates for some of those samples 
cannot be theoretically predicted (and therefore include in a statistical analysis) without a detailed 
simulation of the detector which can only be made by the experimental collaboration itself. Thus our 
results represent the most up-to-date analysis of the atmospheric neutrino data which can be performed 
outside the collaboration. For details on our simulation of the data samples and the statistical analysis 
see the Appendix of Ref . 4l|] . 

For what concerns LBL accelerator experiments, we combine the results on ly^ disappearance from 
K2K 37[ with those obtained by MINOS at a baseline of 735 km after a two-year exposure to the 
Fermilab NuMI beam, corresponding to a total of 3.36 x 10^*^ protons on target 38|]. We also include 
the recent results on — )■ Ue transitions based on an exposure of 7 x 10^*^ protons on target 5ll. l52|. 

In order to test the description of the present data in the absence of 0i3-induced effects we show in 
Fig. [2] the present determination of the leading parameters Arrig^ and 623 for ^13 = and Am|]^ = 
7.6 X 10~^ eV^ from the two atmospheric neutrino analyses and the LBL accelerator results. For 
concreteness we plot only normal ordering; the case of inverted ordering gives practically identical 
results as long as ^13 = 0. This figure illustrates how the bounds on the oscillation parameters 623 
and Aml^ emerges from a complementarity of atmospheric and accelerator neutrino data: |Am|^| is 
determined by the spectral data from MINOS, whereas the mixing angle 623 is still largely dominated 
by atmospheric data from Super-Kamiokande with a best-fit point close to maximal mixing. It is 
interesting to note that there is a very good agreement in the location of the best-fit points from 
SK(I-I-II) and MINOS. This is not the case for SK(I-MI-MII) for which the best-fit point in |Am|i| is 
now lower than the one obtained from LBL. 



2.1.3. Status of 6*13 from global data in 2010 



The third mixing angle 613 is of crucial importance for future oscillations experiments. Fig. [3] summa- 
rizes the information on ^13 from present data, which emerges from an interplay of different data sets. 
An important contribution to the bound comes, of course, from the CHOOZ reactor experiment 54l | 
combined with the determination of lAmg^l from atmospheric and long-baseline experiments. Us- 
ing this set of data, a possible hint for a non-zero ^13 from atmospheric data has been found in 



Refs. 



551 . 156|l . The origin of such a hint has been investigated in some detail in Ref. [481 ] . and more 
recently in l57|. From these results one may conclude that the statistical relevance of the hint for 
non-zero ^13 from atmospheric data depends strongly on the details of the rate calculations and of 
the analysis. Furthermore, the origin of that effect might be traced back to a small excess (at 
the la level) in the multi-GeV e-like data sample in SK-I, which however is no longer present in the 
combined SK(I-|-II) data and is extremely weak in SK(I-|-II-|-III) data. A very recent analysis (ne- 
glecting subleading Am^i effects) from the Super-Kamiokande collaboration finds no evidence of such 



a hint 



34|. 
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Another fragile indication of non-zero ^13 arises from the results of the MINOS experiment. In 
Ref. [39(] the first results on the search for z^^ — )■ v^, transitions were reported, based on an exposure of 
3.14 X 10^*^ protons-on-target in the Fermilab NuMI beam. The collaboration observed 35 events in 
the Far Detector with a background of 27 it 5 (stat) it 2 (syst), corresponding to a 1.5cr excess which 
could be explained by a non-zero value of ^13. Recently a new analysis with double statistics (exposure 



of 7 X 10^*^) has been presented jSll . l52t |. The MINOS collaboration reported the observation of 54 
events with an expected background of 49.1 ± 7.0 (stat) it 2.7 (syst), thus reducing the excess above 
background to 0.7cj. 

An important piece of information on ^13 comes from solar and KamLAND data. The relevant 
survival probabilities are given by 

. cos^ ^is fl — sin^ 201 2 ( sin^ 6)) solar, low energies / KamLAND 
cos^ ^13 sin 612 solar, high energies 

where cf) = /^m^iL/AE and {slv? (p) ~ 1/2 for solar neutrinos. Eq. ([6]) implies an anti-correlation of 
sin^ ^13 and sin^ ^12 for KamLAND and low energy solar neutrinos. In contrast, for the high energy 
part of the spectrum, which undergoes the adiabatic MSW conversion inside the sun and which is 
subject to the SNO CC/NC measurement, a positive correlation of sin^ ^13 and sin^ ^12 emerges. As 



discussed already in [58|,|59[, this complementarity leads to a non-trivial constraint on ^13 and it allows 
to understand the hint for a non-zero value of ^13, which helps to reconcile the slightly different best 
fit points for 612 as well as for Am|]^ for solar and KamLAND separately 4ll. l56l. 1591461 1. 



We found that the inclusion of the new solar data, and in particular of the SNO-LETA results 
tends to lower the statistical significance of ^13 7^ while the results from Vf, appearance from MINOS 
increases it. Within the context of the solar model with higher metallicities (GS98) and for the original 



Ga capture cross-section [49|] , we conclude that the significance of ^13 7^ from solar -|-KamL AND data 
is 79% (I.260") which increases to 81% (1.31(t) after inclusion of the atmospheric, CHOOZ and LBL 
data. We also found that using the solar neutrino fluxes required to fit the lower metallicity data 
(AGSS09) and/or the modified (lower) cross-section for neutrino capture in Ga lowers the best fit 
value of ^13 and its statistical significance. So when using the AGSS09 fluxes and the lower Ga cross- 
section the significance of ^13 ^ from solar -|-KamLAND data is 70% (1.05o") and 76% (1.17o") for 
adding atmospheric, CHOOZ and LBL data. 



2.1.4. Tritium beta decay experiments 



The neutrino mass scale is constrained in laboratory experiments searching for its kinematic effects in 



Tritium /3 decay which are sensitive to the so-called effective electron neutrino mass j62|-l64l] 



22 2,22 2,2 2 

^13^12^1 + Ci3Si2"Z2 + ■5i3m3 , 



(7) 
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At present the most precise determination from the Mainz 65| and Troitsk 66|] experiments give no 
indication in favor of m^^ ^ and one sets an upper hmit 



m^^ < 2.2 eV, 



(8) 



at 95% confidence level (CL). A new experimental project, KATRIN 67|, is under construction with 
an estimated sensitivity limit: m^^ ~ 0.2 eV. 



2.1.5. Neutrinoless double-beta decay experiments 

Direct information on neutrino masses can also be obtained from neutrinoless double beta decay 
(Oz^/3/3) searches provided they are Majorana particles. In the absence of other sources of lepton 
number violation in the low energy Lagrangian, the Oi^/3/3 decay amplitude is proportional to the 
effective Majorana mass of Vf,, ruee, 



me. 



E 



(9) 



which, in addition to the masses and mixing parameters that affect the tritium beta decay spectrum, 
depends also on the phases in the leptonic mixing matrix. The strongest bound from Oz//3/3 decay was 
imposed by the Heidelberg-Moscow group [68] 



ruee < 0.26 (0.34) eV at 68% (90%) CL, (10) 

owever, there are large uncertainties 



which holds for a given prediction of the nuclear matrix element. E 
in those predictions which may considerably weaken the bound 
planned with sensitivity of up to niee ~ 0.01 eV j7o| . 



69( 1 . A series of new experiments is 



2.1.6. The impact of cosmological fits 



Neutrino oscillation data provides as unique information on the absolute neutrino mass scale a lower 
bound 



E 



m, 



> 



y'lAm^ for Normal hierarchy, 



2y|Am3]J for Inverted hierarchy. 



(11) 



Furthermore, neutrinos, like any other particles, contribute to the total energy density of the Universe. 
Furthermore within what we presently know of their masses, the three Standard Model (SM) neutrinos 
are relativistic through most of the evolution of the Universe and they are very weakly interacting 
which means that they decoupled early in cosmic history. Depending on their exact masses they 
can impact the CMB spectra, in particular by altering the value of the redshift for matter-radiation 
equality. More importantly, their free streaming suppresses the growth of structures on scales smaller 



10 



than the horizon at the time when they become non-relativistic and therefore affects the matter power 



spectrum which is probed from surveys of the LSS distribution (see 71[ for a detailed review of 
cosmological effects of neutrino mass). 

Within their present precision, cosmological observations are sensitive to neutrinos only via their 
contribution to the energy density in our Universe, Q^h'^ (where h is the Hubble constant normalized 
to Hq = 100 km s~^ Mpc~^). Qiyh"^ is related to the total mass in the form of neutrinos 

n^h^ = S^/(94 eV) . (12) 

Therefore cosmological data mostly gives information on the sum of the neutrino masses and has very 



little to say on their mixing structure and on the ordering of the mass states (see Ref. T^] for a recent 
update on the sensitivity of future cosmological observations to the mass ordering.) 



In Ref. 73|] we have studied the information on the absolute value of the neutrino mass which can 
be obtained from the analysis of the cosmological data in oc^^CDM + AA^rci + cosmologies where, 
besides neutrino masses, one allows for non-vanishing curvature, dark energy with equation of state 
with cj 7^ — 1 together with the presence of new particle physics whose effect on the present cosmological 
observations can be parametrized in terms of additional relativistic degrees of freedom. To break the 
degeneracies in these models, at least the information from four different cosmological probes must be 
combined. Thus we have performed analysis including the data from CMB experiments, the present 
day Hubble constant HO, measurement, the high-redshift Type-I SN results and the information from 
large scale LSS surveys. 

In Fig. m we plot the 95% allowed regions (for 2 dof) in the planes {rriu^, T.^) and (mee, Sjy). In the 
figure we also show superimposed the single parameter 95% bounds on from different cosmological 
analysis. The figure illustrates the well-known fact that currently for either mass ordering the results 
from neutrino oscillation experiments imply a lower bound on m,^^. On the contrary rriee is only 
bounded from below for the case of the normal ordering while full cancellation due to the unknown 
Majorana phases is still allowed for the inverted ordering. These results show that, even for the 
most restrictive analysis including LSSPS, part of the allowed ranges for rriu^ in the context of the 
owCDM + AA^rei + "T-i/ cosmologies are within the reach of the KATRIN experiment. This is not the 
case for ACDM + models unless only the information of CMB and BAO (or SN) is included. We 
also find that near future neutrinoless double beta decay can test some of the allowed ranges in all 
these scenarios. This will be complementary to the improvement on the expected sensitivity from 



upcoming cosmological probes such as the Planck mission 74i |. 
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Figure 3. Ax^ dependence on sin^ 613 from various data sets as labeled in the figure. The right panel 
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shows Icr ranges AGSS09 and GS98 refer to low and high mctallicity solar models, respectively [53 
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Figure 4. 95% allowed regions (for 2 dof) in the planes (m^^, E^) and {rriee, S^) from the global 
analysis of oscillation data (full regions). We also show superimposed the 95% upper bounds on 
from cosmological constraints for the different analysis as labeled in the figure. 
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2.2. Ois: phenomenology, present status and prospect |l2| 



The leptonic mixing angle ^13 is currently a high-priority topic in the field of neutrino physics, with 
five experiments under way, searching for neutrino oscillations induced by this angle: the reactor 
neutrino experiments Daya Bay 75|, Double Chooz |76l |. RENO 77[ and the accelerator experiments 
NOi^A [78( 1 and T2K [79|. The results of these experiments will be essential for the planning towards 



a possible next generation of long-baseline neutrino experiments able to address leptonic CP violation 
and the neutrino mass hierarchy. 

On the theoretical side, the determination of 613 will provide important information on the mecha- 
nism of neutrino mass generation and the flavour structure in the lepton sector. Considering neutrino 
mass models without any flavour structure, so-called anarchical models, one does expect a value of 
9i3 close to the present bound 80|. If on the contrary experiments would indicate a very tiny value 
for one might wish to have a symmetry reason as an explanation. For example, rather symmetric 
patterns for the mixing matrix are the tri-bimaximal [81[ or the bimaximal 84] mixing matrices. 

The present situation obtained in the global fit of all relevant oscillation data, can be summarized 
according to the updated analysisl of 47|]. The following bounds at 90% (3) CL are obtained: 



( 0.053(0.078) solar + KamLAND 
0.033(0.058) CHOOZ + atm + K2K + MINOS (13) 
0.031(0.047) global data 

The hint for ^13 > coming from the different data sets can be quantified by considering the Ax^ 
for 6*13 = 0: 



sin^ 613 < i 



2.2(1.5) solar + KamLAND 

0.8(0.9) CHOOZ + atm + K2K + MINOS , 

(14) 

0.6(0.7) MINOS i/e appearance 
1.8(1.3) global data 

In tableUare compared the best-fit values for sin^ ^13 and the significance of the hint for ^13 > from 
the global fits to neutrino oscillation data from three different groups. All groups find a non-zero best- 
fit point in the range sin^ ^13 = 0.01 — 0.02. While it is premature to draw strong conclusions from these 
results, upcoming experiments will answer very soon the question whether ^13 is indeed in the range 
indicated by present global analyses. Reactor experiments see a large signal of v'e events, and search 
for a small deviation from the non-oscillation prediction due to 0i3-induced v'e disappearance. These 
are a precision experiment, whose success relies on statistical as well as systematical errors below 
the percent level. Table [H] summarises a few key parameters of reactor experiments. Accelerator 
experiments look for the appearance of the Ve flavour in an almost pure beam, due to oscillations. 



The T2K (TokaitoKamioka) experiment [79(] will use a high intensity off-axis (2.5°) neutrino beam. 



with a peak energy of 700 MeV, generated by a 30 GeV proton beam at J-PARC (Japan Proton 
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reference 



best-fit and la errors significance 



Fogli ct al. [57] 


sin^ 


^13 


= 0.02 ±0.01 


2cr 


Gonzalez-Garcia et al. [Ij (GS98) 


sin^ 6 


'l3 


= 0.0095t°:SJ? 


I.Bct 


Gonzalez-Garcia et al. [Ij (AGSS09) 


sin^ 


^13 


- 0.008i-.^- 


l.lcr 


Schwetz et al. [47] (GS98) 


sin^ 


^13 


= o.oi3tg:°}i^ 


1.5(7 


Schwetz et al. [47] (AGSS09) 


sm 


^13 


= o.oiot-.l^J^ 


1.3(7 



Table I. Comparison of the best-fit values for sin 613 and the significance of the hint for 6*13 > from 



471 include 7 x lO^'^ pot 



different global fits to neutrino oscillation data. The numbers from [l| and 
Ue appearance data from MINOS, whereas [57| is based on 3.14 x 10^° pot. AGSS09 and GS98 refer 
to low and high mctallicity solar models, respectively 53 1. 



Setup 


PTh [GW] 


L [m] 


'T^Dot [t] Events/year Backgrounds/day 


Daya Bay [75] 


17.4 


1700 


80 10 • 10"^ 0.4 


Double Ghooz [76] 


8.6 


1050 


8.3 1.5-10'^ 3.6 


RENO [77] 


16.4 


1400 


15.4 3 • lO'^ 2.6 



Table II. Summary of experimental key parameters of upcoming reactor neutrino experiments. We give 
the thermal reactor power, the approximate distance between reactors and far detector, and detector 
mass, neutrino events per year, and background events per day. all for the far detector. RENO 
backgrounds are the sum of correlated backgrounds as computed in [ttI and uncorrclatcd backgrounds 
as estimated in 1831 . 



Accelerator Research Complex) fired to the Super Kamiokande detector, located 295 km from the 
proton beam target. 

The NOz^A experiment 78( will run at an upgraded NuMI neutrino beam expected to deliver 
6.5 X 10^'^ pot/year, corresponding to a beam power of 700 kW, generating a neutrino beam with 
an average energy E^, ~ 2 GeV and a z^e contamination less than 0:5%. The far detector, placed 
at baseline of 810 km, 14 mrad (0.8°) off-axis, will be a totally active tracking liquid scintillator, 
constructed from liquid scintillator contained inside extruded PVC cells. 

Fig. [5] shows the ^13 discovery reach of the five upcoming experiments expected in 2018. It is clear 
from the figure that the discovery potential of the appearance experiments strongly depends on the 
CP-phase as well as on the neutrino mass hierarchy. We observe that the inverted hierarchy gives a 
weaker sensitivity. Hence, in case no appearance signal is found the final ^13 limit will be set by the 
IH. The different shape of the IH curve for NOz^A results from the anti-neutrino running included in 
the NOi^A run plan. As evident from the figure, reactor experiments are neither sensitive to the value 
of 5 nor to the mass hierarchy. 

The sensitivity of the different experiments to ^13 can be discussed using two different performance 
indicators: the ^13 sensitivity limit and the ^13 discovery potential. The ^13 sensitivity limit describes 
the ability of an experiment to constrain ^13 if no signal is seen. It is basically determined by the 
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Figure 5. Discovery potential of the five upcoming experiments in the plane of sin^ 20i3 and 6 expected 
in 2018, see section 5.3 of Ref. Q for our assumptions on exposure. To the right of the curves a non-zero 
value of can be established at 3ct. For the beam experiments we show normal (solid) and inverted 
(dashed) hierarchies, while reactor experiments are independent of the hierarchy. The four lines for 
Daya Bay correspond to different assumptions on the achieved systematic uncertainty, from weakest 
to strongest sensitivity: 0.6% correlated among detector modules at one site, 0.38% correlated, 0.38% 
uncorrclated among modules, 0.18% uncorrclatcd. 



worst case parameter combination which may fake the simulated ^13 = 0. The ^13 sensitivity limit 
time evolution is shown in Fig. [6l We observe that the global sensitivity limit will be dominated by 
reactor experiments. 

In case of no signal, the 613 limit from beam experiments suffers from the marginalization over the 
CP phase and the mass hierarchy. This situation is very different in case of the discovery potential, 
since there a favourable value of 6 can greatly enhance the sensitivity of the appearance experiments. 
The 013 discovery potentials are shown in Fig. [7] as a function of time. For the beam experiments, 
the dependence on the true value of 5 is reflected by the interval between the solid curves for a given 
time (shaded regions). The dashed curves for T2K and NOz^A correspond to a fixed value for the CP 
phase of 5 = The reactor experiments are not affected by the true 6; the various curves for Daya 
Bay again correspond to the different assumptions concerning systematics as described above. The 
comparison of Figs. [7] and [6] shows that suitable values of 5 may significantly improve the discovery 
potential of beams compared to their sensitivity limit. Indeed, T2K may discover ^13 for smaller ^13 
than Daya Bay in a significant fraction of the parameter space, depending on the achieved systematics 
in Daya Bay. The NOi^A band becomes more narrow due to the complementary information from the 

Evolution of sensitivities under this condition have been shown recently in 

0,0. 
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Figure 6. Evolution of the 6*13 sensitivity limit as a function of time (90% CL), i.e., the 90% CL 
limit which will be obtained if the true ^13 is zero. The four curves for Daya Bay correspond to 
different assumptions on the achieved systematic uncertainty, from weakest to strongest sensitivity: 
0.6% correlated among detector modules at one site, 0.38% correlated, 0.38% uncorrelated among 
modules, 0.18% uncorrelated. 



anti-neutrino running, with the clear disadvantage of being somewhat late. 

In figure [7] is also illustrated how the world sensitivity to ^13 could look like under the assumptions 
of the above schedules and that at each point in time a combined analysis of all available data is 
performed. The discovery reach will be set roughly by the optimal sensitivity of T2K, where the 
reactor experiments play an important role in providing sensitivity for the values of 6 unfavourable for 
T2K. This plot nicely illustrates the interplay between reactor and beam experiments and shows that 
the global reach can be enhanced significantly if experiments of both types are available simultaneously 
with comparable sensitivities. 
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Figure 7. Evolution of the 013 discovery potential as a function of time (3 a CL) for NH, showing 
the global sensitivity reach. The bands for the beams and the global reach reflect the (unknown) true 
value of S. For Daya Bay it is assumed a systematical uncertainty of 0.38% correlated among detector 
modules at one site. 
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2.3. Direct Determination of the Solar Neutrino Fluxes from Neutrino Data Isl 



In 1939, Hans Bethe described in an epochal paper 86(] two nuclear fusion mechanisms by which main 
sequence stars like the Sun could produce the energy necessary to power their observed luminosities. 
The two mechanisms have become known as the pp-chain and the CNO-cycle [87[] . In order to precisely 
determine the rates of the different reactions in the two chains, which are responsible for the final 
neutrino fluxes and their energy spectrum, a detailed knowledge of the Sun and its evolution is needed. 
Standard Solar Models (SSM's) |53|, l88l-l93l] describe the properties of the Sun and its evolution after 
entering the main sequence. 

Till recently SSM's have had notable successes in predicting other observations. In particular, quan- 
tities measured by helioseismology such as the radial distributions of sound speed and density [9Ch93|] 
showed good agreement with the predictions of the SSM calculations and provided accurate informa- 
tion on the solar interior. A key element to this agreement is the input value of the abundances of 
heavy elements on the surface of the Sun 9J]. However, recent determinations of these abundances 



point towards substantially lower values than previously expected |95|, |96[. A SSM which incorpo- 
rates such lower metallicities fails at explaining the helioseismological observations [97[, and changes 
in the Sun mod eling (in particular of the less known convective zone) are not able to account for this 



discrepancy [93, l99 1 . 



So far there has not been a successful solution of this puzzle. Thus the situation is that, at present, 
there is no fully consistent SSM. This led to the construction of two different sets of SSM's, one (labeled 
"GS") based on the older solar abundances [SJ] implying high metallicity, and one (labeled "AGS") 



assumiiig lower metallicity as inferred from more recent determinations of the solar abundances 951.196 1. 



In Ref. 



531 ] the solar fluxes corresponding to such two models were detailed, based on updated versions 



of the solar model calculations presented in Ref. 



Alternatively one may attempt to directly determine the solar neutrino fluxes from the solar neutrino 
data itself. In here we summarize the results of the most up-to-day extraction of the solar neutrino 
fluxes directly from the solar neutrino data from Ref in the framework of three-neutrino oscilla- 



tions. The da ta in cluded compris es th e total rates from the radiochemical experiments Chlorine [23l |. 



Gallex/GNO 100(] and SAGE 25l . Il0[ll |. For real-time experiments in the energy range of neutrinos 



we include the 44 data points of the electron scattering (ES) Super-Kamiokande phase I (SK-I) energy- 
zenith spectrum [26(], the 34 data points of the day-night spectrum from SNO-I [27l |. the separate day 
and night rates for neutral current (NC) and ES events and the day-night energy-spectrum for charge 
current (CC) events from SNO-II (a total of 38 data points) [28j], the three rates for CC, ES and NC 



from SNO-III 29[ , and the 6 points of the high-energy spectrum from the 246 live days of Borexino 32 1 . 
Finally, we include the main set of the 192 days of Borexino data [3l|. Besides solar experiments, we 
also include the latest results from the long baseline reactor experiment KamLAND 35|, |4ll], which 
in the framework of three neutrino mixing also yield information on the parameters Am|]^, 612, and 
^13. In addition, we include the information on ^13 obtained after marginalizing over Am|^, 623 and 
6cp the results from the complete SK-I and SK-II atmospheric neutrino data sets (see the Appendix 
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of Ref. [4l|] for full details on our analysis), the CHOOZ reactor experiment [36|], K2K 37|, the latest 
MINOS Ufj_ disappearance data corresponding to a n exp osure of 3.4 x 10^'^ p.o.t. [33], and the first 
MINOS Ufj_ —7- U(, appearance data presented in Ref. [10l[. 

We do a Bayesian analysis in order to produce the posterior porbability distribution for the param- 
eters (Am|i, 6*12, 6*13, /pp, /tbc, /pep, /i3N, /i5o, /"F, /sB, /hep)- In this model independent analysis we 
assume a uniform prior probability over which we impose a set of constraints, such as the lu mino sity 
constraint which relates the number of neutrinos produced with the total Sun luminosity 102t |. as 
well as those needed to ensure consistency in the pp-chain and CNO-cycle, and some relations from 
nuclear physics. For details on the normalization of the fluxes and the nuclear constraints see 0]. An 
important arises from the 

Our results for the analysis with luminosity constraint are displayed in Fig. [SI where we show the 
marginalized one-dimensional probability distributions for the eight solar neutrino fluxes as well as 
the 90% and 99% CL two-dimensional allowed regions. The corresponding ranges at la (and at the 
% CL in square brackets) on the oscillation parameters are: 



Amii = 7.6 ± 0.2 [±0.5] x 10^^ eV^ , 
sin^ 9i2 = 0.33 ± 0.02 [±0.05] , (15) 
sin^ 6^13 = 0.02 ± 0.012 [1°;°^] , 

while for the solar neutrino fluxes are (in units of cm~'^ s~^): 

-,+0.057 r+0.14i ^ ^10 _ r: ne+0.52 r+l.Si ^ n9 



$pp = 5.9io:--^[;--] X 10- , <i>.Be = 5.08^-- [:ri:g] X 10« 

<I>pep = 1.407tn^^[i°:r7]xl0«, ^r.^ = 7.8^1° [^^^o] X 10^ 



(16) 



$i5o = 4.0t [t^J] X 10^ , $i7F < 5.9 [43] x lO"^ , 

cfsB = 5.02tO;J8 [+0-45] X 10^ $hep =1.3±1.0[t3;0]xl0^ 

All these results imply the following share of the energy production between the pp-chain and the 
CNO-cycle 

-^pp-chain _ nQ^+0.005 r+0.011i , . -^CNO _ n ni /I+0.006 r+0.014i /-, 

— o.98b_o U0.014J —? — - L'-0i4„g Uo.OllJ ' \^') 

in perfect agreement with the SSM's which predict Lq-^q/Lq < 1% at the 3a level. 

As seen in Figs. [8| the inclusion of Borexino has a very important impact on the determination of 
the ^Be, pep and CNO fluxes, and indirectly on the pp flux. 

In order to statistically compare our results with the SSM's predictions we perform two diferent 
significance tests. First we do the analisis without asuming gausianity constructing an statistical 
estimator t from the likelihood, where we found that the GS model has a lower t, tcs = 8.5, while 
tAGS = 11.0 which corresponds to P^f = 43% and P^gs = 20%- 

In the second case we use an estimation of the covariance matrix from the posterior probability 
distribution and we do a chi-square test, we found Xqs ~ ^-^(-^gs^ ~ 74%) versus Xags ~ '^•^ 
(PXgs = 50%). 
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Figure 8. Constraints from our global analysis on the solar neutrino fluxes. The curves in the rightmost 
panels show the marginalized one-dimensional probability distributions, before and after the inclusion 
of the Borexino spectral data. The rest of the panels show the 90% and 99% GL two-dimensional 
credibility regions (sec text for details). 



21 



From these results we conclude that, while the fit shows a slightly better agreement with the GS 
model corresponding to higher metallicities, the difference between the two is not statistically signif- 
icant. This is partly due to the lack of precision of present data. But we also notice that, while the 
measurements of SNO and SK favor a lower ®B flux as predicted by the low metallicity models, the 
determination of the ^Be flux in Borexino and the corresponding determination of the pp flux from 
the luminosity constraint show better agreement with the GS predictions. 

Finally in order to check the consistency of our results we have performed the same analysis without 
imposing the luminosity constraint. This allow us to test the relation between the luminosity of the 
Sun as directly measured with the one infered from the determination of the solar fluxes. 

(neutrino-inferred) ^ ^ ^ q_^^j+o.37] _ (18) 

Lq 

Thus at present, the neutrino inferred luminosity perfectly agrees with the one directly determined 
and this agreement is known with a la uncertainty of 15%. 
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3. A new design for the SPL-Frejus Super-Beam 



In this study we consider a graphite target: this choice constitutes a proven technology in existing 
neutrino beams (i.e. T2K and CNGS). We assume a cylindrical shape with r = 1.5 cm, L = 78 cm 
and p = 1.85 g/cm^. A granular Ti target with the same geometry will also be discussed. 

A new horn model inspired by the one used for the MiniBooNE beam, having a large acceptance for 
forward produced pions, has been adopted giving a reduced contamination from wrong-charge pions. 
The generic layout of the horn is shown in Fig. [9] (Left). Taking advantage of the small transversal 



Figure 9. Left: parametrization of the forward-closed horn. Right: distribution of the figure of merit 
A. See the text for the definition of the samples. 

dimensions, the idea of using a battery of four horns in parallel has been proposed. This arrangement 
allows reduced stress on the targets thanks to the lower frequency pulsing (12.5 Hz) which brings the 
average beam power on each target to a level which is currently considered as a viable upper limit for 
solid targets operations (~ 1 MW). We have verified that placing the horns as central as possible (i.e. 
in mutual contact) causes a minor loss of of the order of 1-2% with a mild loss as a function of the 
radial displacement. 

The approach which was followed in the optimization of the forward-closed horn and the decay 
tunnel uses the final sin^ 2^13 sensitivity, as a guiding principle in the ranking of the system. Given 
the well known dependence on the sin^ 2^13 limit on the 5c p phase, we introduced the figure of merit 
of the focusing system A defined as the (^^p-averaged 99 % C.L. sensitivity limit on sin^ 2^13 in units 
of 10"'^. In a first stage the parameters of the forward-closed horn and of the tunnel were sampled 
with uniform probability distributions within large ranges. The maximal length and radius of the 
horn were limited to 2.5 m and 80 cm in order to maintain a compact design suitable for the operation 
of four horns in parallel. Moreover, the inner radius R\ was limited in [1.2, 4] cm, the lower limit 
corresponding to the "integrated target" limit. At this level the target geometry and the current were 
not varied (/ = 300 kA). The obtained distribution for A is shown in the continuous histogram of Fig. [9] 
(Right). A second scan was performed after fixing the horn inner radius at 1.2 cm and restricting the 
ranges of variation of a set of relevant parameters [J]. The distribution of A for this sample is shown by 
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the dashed histogram of Fig. [9j Finally the horn shape was fixed and a further tuning of the tunnel 
length and radius was performed. The dependence of A on the tunnel variables can be reasonably 
fitted with a quadratic function: A = 0.94 + 2.1 • 10-4(L*""[m] - 31.8)2 + 2.4 • 10-2(i?'™[m] - 2.9)^. The 
distributions of A for L*""^ and i?*"'^ in the neighborhood of the minimum, is shown by the dotted 
histogram of Fig. [9) an improvement of 25-30 % is obtained with respect to the initial sample. Since 
the minimum is relatively broad we chose L*^™ = 25 m and R^^^ = 2 m as central values based on 
practical considerations related to the excavation and shielding of large volumes. This compares to 
the previous values of 40 m of length and 2 m of radius. 

We have also observed that an increase in the current (between 300 and 400 kA) tends to sys- 
tematically produce better sensitivity limits. Data are well fitted by a linear function in (/, Ri): 
A = (9.2 - 0.81 • /[lOO kA])/(7.3 - 0.37 • Ri[cm]). The effect of increasing the current, i.e. a stronger 
magnetic field in the vicinity of the target, is physically equivalent to decreasing the minimum horn 
inner radius. In this way, using a constant //i?i(oc B), allows to roughly work at fixed sensitivity. 

The fractions of v^, i'^, fg and v'e with respect to the total are (98.0%, 1.6%, 0.42%, 0.015%) and 
(95.3%, 4.4%, 0.28%, 0.05%) for the positive and negative focusing modes respectively. In positive 
(negative) focusing mode the (i^e) fiuxes are dominated by muon decays: 82% (90%). The c.c. fluxes 
receive instead a large contribution from K 3-body decays (81 % and 75 % in and "-" focusing 
respectively) with decays from the decay chain of "wronge charge" tt at low energy contributing for 
the rest. These fluxes are available on the internet Q]. 

The discovery potential for ^13 ^ and CP violation improves with respect to the previous design. 
The uncertainty on hadro-production has also been addressed, for the graphite target, at the level 
of sensitivities by exploiting the data of the HARP experiment and different models (FLUKA and 
GEANT4-QGSP). More detailed information on the subject can be found in Q]. 

As undergoing studies in the context of the EUROnu design study have shown technical challenges 
for a solution with an integrated horn-target system, we studied the performances of two additional 
configurations assuming: 1) a graphite target separated from the horn (ST); 2) a granular target 
composed of titanium spheres with diameters of C'(mm) (PB); while keeping the target geometry 
unchanged. In both cases we set for the inner radius of the horn (Ri) a value of 3 cm and for the 
current a value of 350 kA. The 1.5 cm gap between the target and the horn is intended to accomodate 
the cooling infrastructure. Thanks to the favorable surface to volume ratio and the possibility to flow 
transversely the coolant within the interstices of the spheres (i.e. a high pressure flow of He gas) the 
granulat target is expected to have a good behaviour even under extreme irradiation conditions. 

The discovery potentials for ^13 7^ and CP violation for these two conflgurations are compared to 
the ones obtained with the former design based on a mercury target (HQ) and to the performance 
with the integrated target (IT) in Fig. [TUl With respect to the IT design the solution with a separated 
monolithic graphite target and increased current gives limits which are only slightly worsened; the PB 
solution gives practically unchanged performance for 5 > tt and some improvement for 5 < n. The 
granular Ti-target in association with the optimized horn represents then possibly the most appealing 
solution in terms of both physics performance and engineering. 
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Figure 10. The ^13 (Left) and CP violation discovery potential (Right) at 3 a. Known parameters 
were included in the fit assuming a prior knowledge with an accuracy of 10% for 612, 623, 5% for 
Am^i and 3% for Amf2 at 1 tr level. The running time is {2i'+8P) years. Curves are calculated with 
GLoBES 3.0.14. The parametrization of the MEMPHYS Water Cherenkov detector is implemented 
in the publicly distributed AEDL file SPL.glb. 

Bet a- Beams 



Beta-beams are one of the two new beam technologies that have been proposed in the last decade 
to produce intense neutrino beams aiming at distant detectors. Many different beta-beams proposals 
have been studied in recent years. The physics performances of most of them are summarized in 
Sect. ??. The rest of the section covers: atmospheric u,, background at ICAL@INO for a high-7 beta- 
; minimal beta-beams that exploit at most existing european infrastructures [7[; 
an update of the physics potential of electron-captur_e beta-beams [sl; and, the physics potential of a 
high-7 beta-beam within the EU LAGUNA project 



4.1. Performances of Beta-Beam setups as of January 2011 [51 



In this short EUR0nu-WP6 internal note, we review the physics potential of several beta-beam setups 
that have been proposed in recent years, comparing their performances in the sensitivity to sin^ 2^13, 
the CP discovery potential and the sensitivity to the neutrino mass hierarchy. Combination of these 
facilities with other facilities (as the proposed synergy of the SPL super-beam and the 7 = 100 beta- 
beam aiming at the Frejus underground laboratory) and/or with atmo sphe ric neutrino oscillation data 
collected at the same detector of the beta-beam setup under study 103t | are not considered in this 
review. 



4.1.1. "Low"-7: 7 = 100 



Belong to this category all the setups that use existing CERN facilities to boost the ions up to the 
desired energy. In particular, the reference setup is the 7 = 100 ^He/^*Ne beam aiming at a 1 Mton 
class water Cerenkov detector located in the Frejus underground laboratory jl04l . ll05| |. with a baseline 
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Setup 


7 


Ions 


Fluxes [10^8] 


Years 


Baseline 


Detector Technology 


CERN-Frejus, 1 


100 


6He 

i^Ne 


2.9 
1.1 


5 
5 


130 km 


440 Kton WC 


CERN-Frejus, 2 


100 


'^He 
i^Ne 


2.9 X 2 
1.1/2 


2 

8 


130 km 


440 Kton WC 


CERN-Frejus, 3 


100 


^He 

i^Ne 


2.9 X 2 
1.1/5 


2 

8 


130 km 


440 Kton WC 


CERN-Canfranc, 4 


100 


8Li 

8B 


2.9 
1.1 


5 
5 


650 km 


440 Kton WC 


CERN-Canfranc, 5 


100 


8Li 

8b 


2.9 X 2 
1.1 X 2 


5 
5 


650 km 


440 Kton WC 


CERN-Canfranc, 6 


100 


8Li 

8b 


2.9 X 5 
1.1 X 5 


5 
5 


650 km 


440 Kton WC 






8Li 


2.9 


3 






CERN-Canfranc, 4 


100 


8b 

6 He 


1.1 
2.9 


5 
2 


650 km 


440 Kton WC 






8Li 


2.9 X 2 


3 






CERN-Canfranc, 5 


100 


8b 

f^He 


1.1 X 2 
2.9 X 2 


5 
2 


650 km 


440 Kton WC 






8Li 


2.9 X 5 


3 






CERN-Canfranc, 6 


100 


8B 

6 He 


1.1 X 5 
2.9 X 5 


5 
2 


650 km 


440 Kton WC 



Table HI. Summary of the characte ristic s of the 7 = 100 beta-beam setups that have been shown in 
the literature (for a review, see Ref. 109|). 



L = 130 Km. Three options have been considered for this setup, depending on the achievable ion 
rates. A second possibility using the same infrastructures at CERN (i.e., the PS and the SPS) is 
to change the stored ions, going from l ow-Q one s such as ^He/^^Ne to high-Q ones, such as ^Li/^B. 
This possibihty was advanced in Refs. jlOSl . Il07l |. Also in this case, several options depending on the 
achievable ion rates have been considered. In addition to that, it is conceivable to store a mixture 
of low- and high-Q ions aiming at t he sa me baseline, in order to use the first- and second-peak of 
the oscillation probability, as in Ref. [lOSl ] . The characteristics of the low-7 setups are summarized in 
Tab-imi 

In Tab. IIVI we compare the performances of the 7 = 100 beta-beam setups defined above in terms 
of three observables: (1) the minimum value of sin^ 2^13 that can be excluded at 3a in case of a null 
result of a given experiment (sensitivity to ^13); (2) the fraction of the (5-parameter space (known as 
the CP-fraction) for which a non-vanishing 6 can be distinguished by (5 = 0, vr at 3cj (CP discovery 
potential), computed foiI§ sin^ 2^13 = 0.1; and, (3) the fraction of the 5-parameter space for which a 
true positive Am^g can be distinguished by a negative Amfg at 3a (sensitivity to the neutrino mass 



This value is approximately consistent with the recent measurement of 6^13 by reactor experiments 110l - lll2 |. 
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Setup 


(sin2 20i3) . 


CP discovery potential 
CP - fraction for sin^ 26*13 = 0.1 


Sensitivity to sign(Ami3) 
CP - fraction for sin^ 26*13 = 0.1 


CERN-Frejus, 1 


5 X 10"* 


70% 


NO 


CERN-Frejus, 2 


6 X lO""* 


70% 


NO 


CERN-Frejus, 3 


1 X 10"^ 


61% 


NO 


CERN-Canfranc, 1 


1.5 X 10-^ 


58% 


51% 


CERN-Canfranc, 2 


7 X 10-'' 


72% 


61% 


CERN-Canfranc, 3 


2 X 10"* 


78% 


100% 


CERN-Canfranc, 4 


1.7 X 10-3 


66% 


100% 


CERN-Canfranc, 5 


7 X 10-4 


71% 


100% 


CERN-Canfranc, 6 


3 X 10-4 


79% 


100% 



Table IV. Summary of the performances of the 7 = 100 beta-beam setups at 3a in terr tis of : sensitivity 
to 6*13; CP discovery potential; sensitivity to the neutrino mass hierarchy. From Ref. jiod ]. 



hierarchy), computed for sin^ 2^13 = 0.1. Notice that in all cases we have considered no atmospheric 
neutrino background. This means that we have always assumed that the duty cycle at which the 
beta-beam is operated is tight enough to make this background negligible. The duty cycle for which 
this approximation is valid differs depending on the setup. How the performance are deteriorated 
when the duty cycle is relaxed is not shown in this table. The impac t of the atmospheric neutrino 
background on the beta-beam performance has been studied in Ref. 1091] for the low-7 case. See 
Sect. 14.21 in this report for an estimate of the effect of atmospheric neutrino background for high-7 
setups. 



4.1.2. "High"-7: 7 > 350 



In Ref. 



113l . Ill4i | , the possibility of using an upg rade of the SP^^ to boost radioactive 



ions up to 



higher 7 values. In particular, it was shown that the sensitivity to ^13 and the CP discovery potential 
of a beta-beam with ^He/^^Ne ions boosted at 7 = 350 was extremely good and competitive with 
the Neutrino Factory in some part of the parameter space. As for the 7 = 100 option, the beam 
was directed towards a 1 Mton class water Cerenkov detector, located this time at the underground 
laboratory of Canfranc (L = 650 Km from CERN). This option was later adapted for the Gran Sasso 
underground laboratory, where such a big detector cannot be hos ted, by r eplacing it with a 100 Kton 
iron detector (that can be magnetized, to reduce backgrounds) 115l4ll7l | or with a 50 Kton TASD 
1181]. It was later noticed that, when using high-Q ions with hig h-7, the n eutrino flux is peaked 



around the resonant energy for — ?• conversion in Earth matter [I19l . Il20t ] . This makes a high-Q 
high-7 beta-beam aiming to a very far 50 Kton iron detector (at L ~ 7000 Km) an extremely good 
experiment to measure the neutrino mass hierarchy. As a consequence, several two-baseline beta-beam 



The SPS+, proposed at that time within the of maintenance and upgrade programme of the LHC. 
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Setup 


7 


Ions 


Fluxes [10^8] 


Years 


Baseline 


Detector Technology 






6 He 


2.9 


5 






High-7, 1 [113] 


OOU 






700 km 


500 Kton WC 




i^Nc 


1 1 
1 . i 













6 He 


2.9 


2 






High-7, 2 [118] 


350 




700 km 


50 Kton TASD 




iNe 


1 1 
1.1 


Q 
O 






Two baselines, 1 [1211 




8Li 


3 


5 


2000 km 


50 Kton MIND 


350 


8b 








3 


5 


7000 km 


50 Kton MIND 


Twn ViaQplinpc; 9 f1911 


oo\j 


8 Li 


5 


5 


2000 km 


50 Kton MIND 










8b 








7nnn 
( uuu Km 


ou ivLon iViiiNu 


1 WO Dtiseiines, o 




8t ; 


1 n 
lU 


r 



zUUU Km 


oU Kton iViliNU 


oOU 


8B 






1 n 
lU 





(Uuu Km 


ou ivLon iViiiNU 




390 


8Li 


0.6 X 3 


2.5 


7000 km 








^B 






50 Kton MIND 


Cocktail, 1 [1221 


656 


0.6 X 3 


2.5 












350 


^He 


3 


2.5 
















650 km 


500 Kton WC 




350 


3 


2.5 








390 


^Li 


3 


2.5 
















7000 km 


50 Kton MIND 




656 


8B 


3 


2.5 






Cocktail, 2 [1241 










575 


6 He 


3 


2.5 










650 km 


50 Kton TASD 




575 




3 


2.5 







Table V. Summary of the characteristics of the high-7 beta-beam setups that have been shown in the 
literature. 



Setup 


('sin2 26ii3) . 


CP discovery potential 
CP — fraction for sin^ 29i3 = 0.1 


Sensitivity to sign(Am^3) 
CP - fraction for sin^ 26*13 = 0.1 


High-7, 1 


1.6 X 10-'' 


93% 


100% 


High-7, 2 


5 X 10-'' 


75% 


85% 


Two-baselines, 1 


1 X 10-3 


44% 


100% 


Two-baselines, 2 


3 X 10-* 


62% 


100% 


Two. baselines, 3 


1.5 X IQ-'' 


74% 


100% 


Cocktail, 1 


1.8 X IQ-'' 


81% 


100% 


Cocktail, 2 


5 X lO-'' 


73% 


100% 



Table VI. Summary of the performances of the 7 > 350 beta-beam setups at 3cr in terms of: sensitivity 
to 013 ; CP discovery potential; sensitivity to the neutrino mass hierarchy. 



setups have been proposed, using a resonant beam aiming to L = 700 Km and a second beam aiming 
to a m oder ate distance det ecto r (L ~ 2000 Km when using Li/B [1211], L = 650 Km when using 
He/Ne |122l |: see also Refs. 123l . Il24l |). These beams was shown to be competitive with the Neutrino 
Factory in most part of the parameter space. The characteristics of the high-7 setups are summarized 
in Tab. |Vl In Tab. IVII we compare the performances of the high-7 beta-beam setups defined above. 
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Figure 11. Atmospheric neutrino events in 50 Kton ICALQINO detector in 5 years. 



4.2. Atmospheric neutrino events at ICAL@INO and high Q /3-beam [61 



1. Atmospheric events: 

• First of all, we compute the expected number of atmospheric events in 50 kton ICAL@INO 
detector in 5 years. We have calculated separately the neutrino and anti-neutrino events. 



l^^i + I'e ^ 



oscillation channels. In case of 
Vu)- Fig. 1 shows the expected 



For neutrinos, we have considered {u^ 
anti-neutrinos, the considered channels are (z^^ — )• z/^ -|- 
number of events as a function of neutrino energy at ICAL@INO detector in 5 years. Here 
we have done the integration over 9 and (p in their entire range. All other details of the 
simulation have been mentioned on the body of the figure itself. Next, we will present these 
number of events in a tabular form (Table. 1 and 2). 

For beam studies, the atmospheric events which will occur along the beam direction and 
arount it will serve as background. The zenith angle for CERN-INO baseline is 124° and in 
ICAL@INO detector, the angular resolution will be around 15° at most and with higher 
energies, the angular resolution improves a lot. Therefore, in our next study, we have 
considered a zenith angle range of 109° to 139°. In this zenith angle range, the atmospheric 
events that you expect at 50 kton ICAL@INO detector with 5 years of data taking is given 
in Table. 3. 

One can see from Table. 3 that the atmospheric neutrino flux falls steeply with energy and 
is expected to produce much fewer events for the energy range that we are interested in for 
CERN-INO beam study. Therefore, we need to see that how sensitivity will be affected in 
CERN-INO /3-beam set-up with the increase in threshold. 
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Energy Bins 


V events 


V events 


GeV 




(pp. Vf_, + i>e ^ D^j) 


0.4 - 0.6 


1474 


418 


0.6 - 0.8 


1251 


370 


0.8 - 1 


989 


314 


1 - 1.2 


807 


273 


1.2 - 1.4 


653 


236 


1.4 - 1.6 


562 


209 


1.6 - 1.8 


479 


185 


1.8 - 2 


409 


160 


2 - 2.2 


354 


139 


2.2 - 2.4 


308 


126 


2.4 - 2.6 


271 


112 


2.6 - 2.8 


237 


102 


2.8 - 3 


212 


89 


3 - 3.2 


189 


80 


3.2 - 3.4 


170 


72 


3.4 - 3.6 


153 


64 


3.6 - 3.8 


140 


61 


3.8 - 4 


127 


54 


4 - 4.2 


114 


51 


4.2 - 4.4 


106 


46 


4.4 - 4.6 


97 


42 


4.6 - 4.8 


89 


39 


4.8 - 5 


83 


36 



Table VII. Atmospheric neutrino events in 0.4 to 5 GeV range at 50 Kton ICALQINO detector in 5 
years. 

Fig. 2 shows that the sensitivity of the experiment remains almost stable against the varia- 
tion of the threshold energy upto 4 GeV. It means that we can work with a threshold of 4 
GeV or so and in that way, we can reduce the atmospheric background a lot as can be seen 
from Table. 3. 

• The fact that INO has charge identification capability further reduces the atmospheric back- 
ground. The most important handle on the reduction of this background comes from the 
timing information of the ion bunches inside the storage ring. For 5T magnetic field and 
7 = 650 for '^B ions, the total length of the storage ring turns out to be 19564 m. We have 
checked that with eight bunches inside this ring at any given time, a bunch size of about 
40 ns would give an atmospheric background to signal ratio of about 10~^, even for a very 
low sin^ 2^13 of 10^'^. For a smaller bunch span, this will go down even further. In addition, 
atmospheric neutrinos will be measured in INO during deadtime and this can also be used 
to subtract them out. 
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Energy Bins 




V events 




V events 


GeV 




-> I'/j + J^e - 










5 - 5.2 






78 






34 


5.2 - 5.4 






70 






31 


5.4 - 5.6 






67 






29 


5.6 - 5.8 






63 






28 


5.8 - 6 






58 






26 


6 - 6.2 






54 






24 


6.2 - 6.4 






51 






22 


6.4 - 6.6 






48 






21 


6.6 - 6.8 






45 






20 


6.8 - 7 






42 






19 


7 - 7.2 






41 






18 


7.2 - 7.4 






38 






17 


7.4 - 7.6 






36 






16 


7.6 - 7.8 






34 






15 


7.8 - 8 






33 






14 


8 - 8.2 






31 






14 


8.2 - 8.4 






29 






13 


8.4 - 8.6 






29 






13 


8.6 - 8.8 






27 






12 


8.8 - 9 






26 






12 


9 - 9.2 






25 






11 


9.2 - 9.4 






24 






11 


9.4 - 9.6 






23 






10 


9.6 - 9.8 






22 






10 


9.8 - 10 






21 






9 



Table VIII. Atmospheric neutrino events in 5 to 10 GcV range at 50 Kton ICAL@INO detector in 5 
years. 



Energy range (GeV) 


Total V events 


Total V events 


1 - 12 


214 


94 


2 - 12 


155 


69 


3 - 12 


114 


48 


4 - 12 


90 


39 



Table IX. Atmospheric neutrino events at 50 Kton ICAL@INO detector in 5 years in the zenith angle 
range of 109° to 139°. Here full integration has been done over 0. 



4.3. A minimal Beta Beam with high-Q ions to address leptonic CP violation |0] 



The Beta-Beam concept and its different energy configurations have been discussed in details in 



104 



1251] and References therein. Here we focus on a Beta Beam designed with the aim of leveraging at 
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Figure 12. It shows the effect of changing the threshold energy of the detector. The sensitivity of the 
experiment remains almost stable against the variation of the threshold energy upto 4 GeV. Here we 
have used and ^Li ions with 5 years of running and we h ave assumed 1.1 x 10^® (2.9 x 10^^) useful 
ions decays per year. This figure has been taken from jl2d |. 



most existing facilities and, in particular, the CERN acceleration complex. Beta Beams that are able 
to accelerate radioactive ions to high energies and produce multi-GeV z/g and z^e allow for the use of 
high-density detectors, which, in turn, might be hosted in moderate-size underground laboratories. 
For a CERN-based Beta-Beam, the natural option to host the far detector is a laboratory located 
at a distance 0(600 — 700) Km from the neutrino source. The facility that exploits at most existing 
European infrastructures is a multi-GeV Beta-Beam based on the CERN-SPS accelerator pointing to a 
massive, high-density detector located in one of the experimental halls of the Gran Sasso laboratories 
0]. The next cheapest alternative could be represented by the Canfranc Underground Laboratories 
in Spain, where some engineering would be however needed (albeit not so impressive as for a Mton 
class Water Cerenkov detector). 



Beside the huge practical interest o f ex 
additional infrastructure investment 



[ploiting in an optimal manner all European facilities without 
126t |. this detailed assessment is particularly relevant at present 
times: since 2009, machine studies for the Beta-Beam are concentrated on facilities that accelerate 
ions with Q-values larger than originally envisioned {Q ^ 13 MeV for ^Li and ^B, to be compared 
with Q 3 MeV for the ions considered in the original design, ^He and ^^Ne) using the existing 
SPS machine 



2C 



127l |. This option ideally fits the "minimal" scheme mentioned above provided 



that neutrinos are pointed toward the underground halls of LNGS. Other options either based on 
low density detectors and/or on new terminal boosters at larger energies tha n the CERN-SPS have 



also been studied in literature: for details, we refer the reader to Ref. 
Refs 



105|,|l2a,ll29|] for low-Q ions accelerated by the SPS, Refs. [113l . lll4l . lll6l4ll8l . 



125 ] and, in 



Beta Beams (using facilities different from the SPS to accelerate ions) and Refs. 



13C 



loe 



particular, to 



131 



108 



for h igh 



lOS, 



118- 



32 



124 



132 



135l ] for high-Q Beta Beams (either at low-7 and high-7). 



The facihty that we consider here does not differ from the basehne EUROi^ [l09l . 11361 ] design, but 
for a high density far detector located in a pre-exist ing h aU at LNGS (for details on the facility and 
on the detector set-up and simulation we refer to 0, Ill6| . 

The physics performance of the proposed setup have been studied in terms of two observables, 
defined as follows: 

the CP discovery potential: for a given point in the parameter space, we will say that CP violation 
can be discovered if we can rule out the no CP violation hypothesis {6 = 0° and 180°) at 3a 1 
d.o.f., after marginalizing over all the remaining parameters for both possible hierarchies, 
the sgn(Am23) reach: this is defined as the region of the (sin^ 2^13, (5) plane for which the wrong 
hierarchy can be eliminated at 3a. Below this value of sin^(20i3), the predictions for the wrong 
hierarchy cannot be distinguished from the data corresponding to the right hierarchy, at a 
statistical significance of 3a. 
For a far detector of 100 kton mass, a /3^-emitters (^B) flux of approximately 6 x 10^^ useful decays 
per yeai0 is needed to observe CP violation in a large fraction of the parameter space (60%) for any 
value of 9is that gives a positive signal at T2K (^13 ^3°). This sensitivity to 6 is deteriorated for 
(5 < due to the occurrence of the vr-transit, as observed in other facilities. The '^B flux must be 
accompanied by a ^Li fiux of ~ 3 x 10^^ decays per year. Present studies on the ionization cooling 
technique or on ISOL-type targets indicate that such a large ^Li fiux could be feasible. Moreover, 
the former technique should produce f3~^ and (3~ emitters at a similar rate although ^B ions interact 
stronger than ^Li ions with materials in the target and in the recollection region. To achieve the fiuxes 
above clearly represents the most challenging task for accelerator R&D but it is a viable option with 
respect to ^^Ne, where ISOL-type targets fall almost two orders of magnitude short of the goal. 

In the same configuration, we find a non-negligible sensitivity to the neutrino mass hierarchy that 
extends up to ^13 ~ 4° for positive (negative) values of 5 for normal (inverted) hierarchy. In the 
opposite parameter area, i.e. for negative (positive) values of 6 and inverted (normal) hierarchy, the 
combination with atmospheric data collected during the Beta Beam run by the same magnetized 
detector further improves such sensitivity at large ^13 (~ 6°). Combination of atmospheric data 
with Beta Beam-driven ones should also be able to solve part of the 7r-transit deterioration discussed 
above. 



4.4. Update on the physics of Electron Capture neutrino beams 



Electron capture is a decay channel available to proton-rich nuclei that competes with positron decay 
in a manner that depends on the decay Q-value. Such decay channels are sources of mono-enegetic 
neutrin os s i nce there is no electron in the final state with which to share energy. It has been sug- 
gested 



137 



140l | that ions with dominant electron capture channels can be used as an alternative to 



the standard Beta Beam ions to instead source an intense and collimated beam of mono-energetic 
* This is about three times the flux proposed for ^**Ne, where Fo ~ 2 x 10^*. 
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Figure 13. (Left panel) A discrete approximation to the integrated CP-fraction for a mono-energetic 
neutrino beam, sourced at CERN and directed at a 440 kton Water Cerenkov detector located at the 
Canfranc laboratory, as a function of the choice of boost pairings chosen. The plot shows the results 
for a nm time ratio of low boost : high boost = 7:3. (Right panel) CP- violation discovery plots for 
the boost choices (71,72) — (450, 150) (red solid lines) and (280, 160) (blue dashed lines). 



neutrinos. A neutrino beam of this type could be used to accurately probe the energy dependence of 
the z^e T^fi appearance probability: a means to observe CP-violation without the need for an anti- 
neutrino channel. However, since a single neutrino energy is insufficient to determine the unknown 
neutrino mixing parameters and resolve any degeneracy, it is necessary to includ e at least two ion 
boosts in any physics strategy. In a recent study 0], the previous works 137 . 138 1 were extended to 
examine the physics reach of a ^^''Dy electron capture beam sourced from CERN and directed towards 
a 440 kton fiducial mass Water Cerenkov detector 650 km distant to include the study of the choice 
of boost, the relative run times of each boost, the number of useful ion decays, and the number of 
atmospheric neutrino events. 

In the ffi'st instance, the number of useful ion decays was fixed to 10^*^ per year, and atmospheric 
neutrino backgrounds were neglected so that the effect of altering the ion boosts and run times could 
be investigated. The results of this initial analysis is presented in the left panel of Fig. [T3l for a choice 
of ion boosts (71,72), the percentage of the parameter space sin^ 2^13 G (10~^, 10~^), 5 G (—180, 180) 
for which the existence of CP-violation can be demonstrated (referred to hereafter as CP-coverage) at 
99 % confidence level is shown. (The procedure is outlined in detail in Q].) The results presented are 
for a low boost run taking 70 % of the whole experimental run. All results in this section use 2 degrees 
of freedom statistics. Little variation was found for the relative run times; however, two regions of the 
(71,72) plane returned large CP-coverage. The naive pairing of (almost) first and second oscillation 
maximum had a CP-coverage of approximately 35 % but was asymmetric in 5 — ?• —5 (right panel, red 
solid lines). For 5 < 0, this choice could rule out CP-conservation down to sin^ 2^13 = 10~^, but only 
sin^ 2013 ~ 2 • 10"'^ for 5 > 0. There exists a slightly larger coverage for the pair (280, 160). This 



34 




Figure 14. CP- violation coverage for (71,72) = (280,160) (left) and (71,72) = (440,150) (right) as a 
function of the number of useful decays and number of atmospheric backgrounds per year; both with 
71 for 70 % of the experimental run time. 



choice returned a (roughly) symmetrical CP-sensitivity region with CP-conservation ruled out down 
to sin^ 2^13 ~ 4 • lO"'* (right panel, blue dashed lines). 

However, no R&D has been carried out on the feasibility of sourcing 10^^ us eful decays per year for 
high proton number ions. A brief study adapting EURISOL Beta Beam codes 14ll | indicated that the 
decays are too slow and large tune-shifts would need to be accommodated. With existing technology, 
the number of useful decays is two orders of magnitude too small. At best, therefore, 10^^ useful decays 
per year appears to be overly optimistic and should be considered as the hard limit on the yearly rate. 
To take this experimental fact into account, the ability to uncover CP-violation was investigated, for 
the two boost pairings identified above, as a function of the number of useful ion decays and the 
level of atmospheric neutrino background. The results are summarised in Fig. [T^ again using the 
CP-coverage as an indicator. In both cases, the results were more volatile to small changes in the 
useful decay rate than the atmospheric background. It was found that drops in the total event rate 
can lead to degenerate solutions dropping below the required statistical significance of the test thus 
reducing the sensitivity in a manner beyond simply a scaling of the 01- The background merely 
reduces the without interfering with its ability to rule out degenerate solutions. Its manifestation 
is simply to push the sensitivity contours inwards. 

The feasibility of an electron capture beam remains open question; however, the most optimistic 
parameterisation considered in Q], and shown in Fig. [TBI is likely to be beyond the hard limits imposed 
by available technology. It must be conceded that if ^13 were to be very small, the electron capture 
will not be a competitive facility, unless a fast decaying ion that can be produced in large quantities 
is found. This can be seen explicitly in Fig. 1151 where the CP-discovery for the boost pairs (280,160) 
and (440,150), and 20 atmospheric events per year, are compared to the 7 = 350,350 Beta Beam 
introduced in [114f | with the baseline shortened to CERN-Canfranc, as in this study. The simulation 



used 1.1 X 10^^ useful ^^Ne decays per year and 2.9 x 10^*^ useful ^He decays year with an equal split 
between neutrino and anti-neutrino running for the Beta Beam. 
We see that for a (440,150) electron capture facility, the physics return is poorer but relatively com- 
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Figure 15. CP-violation sensitivity at 99 % confidence level for iVatm — 20 per year for (71,72) = 
(440,150) (black solid) and (280,160) (blue dot-dashed). The red dashed line displays the physics 
return for a standard Beta Beam directed along the CERN-Canfranc baseline with a boost 7 = 350 
for both neutrinos and anti-neutrinos. 



petitive for b < 0°. For these S, the (280,160) facility is only competitive for sin^ 2^13 G [10^^, 10"^]. 
For S > 0°, the (440,150) facility is outperformed by the standard Beta Beam by 2 orders of magnitude 
whilst, again the (280,160) facility is only competitive for sin^ 2^13 G [10~^, 10^^]. The physics reach 
of the electron capture machine is ultimately limited by the need to run for substantial periods at low 
energies where the appearance event rate is small. The availability of a mono-energetic anti-neutrino 
beam could resolve this issue. Ions that decay through bound beta decay could source such a beam; 
however, this is not possible from a practical point of view with available technology [s, 142]. 



4.5. High-7 Beta Beams within the LAGUNA design study j3| 



Within the LAGUNA design study 143l |. seven candidate sites are being assessed for their feasibility 
to host a next-generation, very large neutrino observatory. Such a detector will be expected to feature 
within a future European accelerator neutrino programme: Superbeam, Beta Beam or, if the detector 
is magnetised, a low energy Neutrino Factory. In (9|, a high boos t ^^N e and ^He Beta Beam sourced 



at CERN and directed towards a 50 kton Liquid Argon detector [144i | located at the LAGUNA sites: 



Slanic (L = 1570 km) and Pyhasalmi (L = 2300 km) was investigated. The study assumed the 



availability of an 1 TeV SPS, as envisaged in some LHC upgrade options 145l |. Using boosts of 
7 = 350 for both neutrinos and anti-neutrinos run for 5 years each, the ability to distinguish ^13 from 
zero, rule out CP-conservation, and determine the correct neutrino mass ordering was tested for both 
baselines. The simulations were then repeated but with the boost for ^^Ne increased to 7 = 570. At 
present, no detector response data is available for a Beta Beam neutrino flux incident on a large Liquid 
Argon detector; however, they are expected to possess excellent energy resolution for appearance 
events and background reduction capabilities whist maintaining a high detection efficiency. The results 
presented here take this optimism at face value and therefore should be considered near the limit of 
the possible physics reach of the two setups. 
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Figure 16. 3cr C.L. contours for discovery for non-zero ^13 (left), CP-violation (centre), and sign(Am§2) 
determination (right). In each plot, the solid line corresponds to a 7 = 350,350 ^^Ne and ^He Beta 
Beam for a 50 kton Liquid Argon detector located at Slanic {L = 1570) (top line); and Pyhasalmi 
{L = 2300 km) (bottom line). The dashed lines correspond to a 7 = 570,350 ^®Ne and ^He Beta 
Beam. 



The physics reach of the two baselines considered in summarised is Fig. [1^1 The results presented 
here have been simulated assuming 3 x 10^^ useful decays per year for each ion. An energy threshold 
of 0.4 GeV is taken with energies up to 4 GeV considered with 12 energy bins and a dE/E = 15 %. 
f^-appearance detection efficiency is flat at 80 % and neutral current backgrounds are taken at 0.1 % 
of the unoscillated flux. The LAGUNA study assumes a Liquid Argon detector of mass up to 100 
kton. 50 kton is taken here to bring the simulations in line with other studies in the literature and 
to introduce some redundancy in the event that the parametrisation used here is too optimistic (e.g 
1.5 X 10^^ useful decays per year instead of 3 x 10^^). 

The best sensitivity to non-zero ^13 and CP-violation is found for the CERN-Slanic baseline. This is 
to be expected since, with the same source, the flux is larger for this baseline. The weaker matter effect 
means that the sign(A?n3]^) degeneracy interferes less with these measurements. For the 7 = 570, 350 
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boost pair, non-zero ^13 can be seen down to siv?26 ~ 10~^, and, for both boost pairs, at all values 
of 5 for sin^ 2^13 > 10~^. However, there is a marked difference between the two boost pairings for 
sensitivity to CP- violation. For the 7 = 350, 350 pair, the ability to rule out (5 = 0° or 180° is restricted 
to sin^ 2^13 > 10~^, but increasing the boost of the ^®Ne ions returns a large region of parameter space 
for (5 < 0° and centred on sin^ 2^13 = 10~^. There is little enhancement on the region for low boost 
pairing. This is suggestive that degeneracy is a problem for the 7 = 350, 350 boost pairing; especially 
for sin 26*13 ~ 5 • 10-3 - 1 • 10-2 where there is a gap in CP-violation sensitivity. The lower event rates 
for the longer baselines mean that the data is insufficient to reduce the significance of some degenerate 
solutions. Although, the ability to rule out the incorrect mass ordering is poor relative to Neutrino 
Factories, for a Beta Beam it is not intrinsically bad. For the high boost run, the correct ordering can 
be indentified down to sin^ 2^13 = 2-10-'^, with determination for all values of 5 for sin^ 2^13 > 10"^ in 
both cases. Although the increase in the ^'^Ne boost improves the reach, it does not do so significantly. 
Increasing the boost makes data from higher enegies available without significantly altering the event 
rate and composition at lower energies. Since European baselines make use of low and high energies 
in combination, improving the event rate in one region without the improving other need not, and has 
not, dramatically improved the physics return. The low event rate at low energies is still insufficient 
to break the degeneracy for small values of 6*13. 

The physics return for the CERN-Pyhasalmi baseline is weaker for each of the physics indicators, 
with little sensitivity for sin^ 2^13 < 10-^. Principally, this is due to the L~'^ dependence of the un- 
oscillated neutrino flux. In particular, the ability to determine the correct mass ordering is very poor 
even given the large matter effect at this baseline. The true and incorrect mass ordering solutions will 
be sufficiently separated in (^13, (5) space; however, the low event rate means that the solution regions 
at 3fT will be large and possibly merged together. When combined with a large solution region from 
the low energy bins (also owing to low event rates), the data is insufficient to break the degeneracy 
for small ^13. 
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5. Physics at the Neutrino Factory 



Neutrino factories have been studied as a mean to produce intense and collimated neutrino beams 
aiming at far detector(s) for the last ten years. Their impressive physics potential make them as 
the ultimate neutrino oscillation facility to improve our precision on known leptonic mixing matrix 
parameters and to measure the few remaining unknowns. The main activity in the field in 2010 focused 
on the possibility to "stage" a Neutrino Factory, so that at each upgrade a new physics program me 
is at hand [lO|; on the physics potential of a "low-energy" Neutrino Factory that could optimize its 



capability to measure leptonic CP violation in case 613 were large (as we now know it is the case); 
and, the study of an important source of background to the "golden channel" — )• oscillation that 
was previously overlooked (the r-contamination problem, see [l^). 



5.1. Neutrino Factory in stages |lO 



Here we report on neutrino oscillation physics with a neutrino factory in stages, Ref. We include 
the possibility of upgrading the muon energy, adding another baseline, and increasing the detector 
mass within the same program. Two aspects are taken into account depending on the possible ^13 
discovery by the next generation of experiments. For the large ^13 case (sin^ 2^13 > 10~^), on the on e 



hand, a low energy neutrino factory (LENF) with ~ 4GeV to 5 GeV, discussed in Refs. [146l4l49l |. 
could be sufficient to discover CP violation (CPV) and the mass hierarchy (MH). We perform an 
optimization study of the physics reach, and we identify the "minimal" neutrino factory. On the other 
hand, we propose a staging scenario to chase ^13 if ^13 is smaller, i.e., sin^ 2^13 < 10"^. For details, 
see Ref. Q. 

Minimal neutrino factory fo r th e large ^13. We investigate the minimal requirements for a 



neutrino factory suggested in Ref. [13J]: 

1. 5a independent confirmation of sin^ 2^13 > (for any 6cp)- 

2. 3(7 determination of MH for any (true) Sqp- 

3. 3cr establishment of CPV for a certain fraction (such as 80%) of all (true) 6cp- 

We determine the minimal necessary luminosity as a function of baseline with respect to the above 
performance indicators in Fig. [T71 where the requirement for ^13 can typically be easily met. The 
luminosity is expressed by a scale factor (SF) where SF=1 corresponds to 2.5 • 10^*^ useful muon 
decays per year, which is the IDS-NF standard per baseline and polarity. For the LENF, the minimal 
baseline is determined by the MH reach, and the minimal SF by the CPV reach. The nominal 
luminosity SF=2.8, coming from the injection of all muons in the same storage ring and the re- 
optimization of the frontend, is sufficient for the CPV measurement for 80% of all true 5cp and for 
the MH measurement for all 6cp in the baseline window 1100 km ^ L < 1400 km for both values of 
sin^ 2^13 (left and right panel) . One can read off these figures that luminosity is clearly an issue for 
large sin^ 2^13. If, for instance, only a lower SF can be achieved, the CPV discovery reach decreases 
accordingly. The "minimal" (optimal) LENF, i.e., the one with the lowest SF, is in both panels at 
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Figure 17. Discovery of CPV (dark/red) and MH (medium gray/light blue) for the one baseline 
(minimal) LENF as a function of baseline and luminosity scale factor SF. Discovery reach is given 
within the shaded regions at the Scr CL. where for CPV a fraction of Sep of 75% or 80% is required (as 
indicated), and for the MH a fraction of 5cp of 100%. The stars show the baseline with the minimal 
SF: in the left panel (1100 km, 2.0) and in the right panel (1150 km, 2.6). The nominal luminosity of 
the LENF is given by SF=2.8. Here the true value of sin^ 2^13 is chosen as given in the plot panels, 
and a normal hierarchy is assumed. The matter density uncertainty is assumed to be 2%. Figure taken 
from Ref. 0. 



about L ~ 1100 km, which is close to the FNAL-Homestake baseline L = 1290 km. 

Staging for the small ^13. After several years of data taki ng fr om Daya Bay and the other next 
generation experiments, we will know whether sin^ 2^13 < 0.01 150l |. What has to be done if ^13 has 



not been found by these experiments? Here we propose a plausible staging scenario, starting with 
a low energy neutrino factory. Meanwhile, it is open to take three upgrade options increasing the 
muon energy, adding another baseline and increasing the detector mass. We proceed in three phases 
of data taking, five years each. As we show in Fig. [T2] (5ct), phase I represents a low energy neutrino 
factory, phase II includes the energy upgrade to a high energy neutrino factory, and phase III considers 
additional upgrades, such as a larger detector or an additional baseline. In any phase, A combination 
of the data with the preceding phase is considered by default. We see the neutrino factory in stages 
for small sin^ 2^13, where a second baseline is used in phase III (solid curves) or the detector upgrade 
(dashed curves). At the 5a confidence level, the strength of the magic baseline to resolve degeneracies 
becomes more important, especially for the mass hierarchy discovery reach. The largest increase in 
the discovery reaches just comes from the energy upgrade. 

In summary, a realistic program may include components of LENF and HENF, even if sin^ 2^13 is 
small. 
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Figure 18. 6*13, MH, and CP discovery reaches for a neutrino factory in three phases: Phase I (hght 
gray/green) is a low energy neutrino factory with a magnetized TASD, phase II (medium gray/red) 
adds an energy upgrade with a MIND at the 4000 km basehne, and phase III (dark gray/blue) includes 
another (magic) baseline (solid curves) or a detector upgrade at the 4000 km baseline (dashed curves). 
5(7 confidence level. Figure taken from Rcf. [lo| . 



5.2. LENF Overview [111 



The low-energy neutrino factory (LEN F), which uses muons of energy ~ 5 GeV and a baseline of 



1000km, was first proposed in Ref. 



146 



143]. Since then, developments to the accelerator and 
detector designs have enabled the experimental simulations to be refined and detailed optimisation 
studies to be performed [ll[. The key finding was that given sufficiently high statistics, an optimised 
LENF can have excellent sensitivity to the standard oscillation parameters, competitive with, and 
even better than, the high energy neutrino factory for large ^13 (sin^ 2^13 > 10^'^). The possibility of 



observing the p latin um channels (z^^ — )• and u^j 



and Vp 



in addition to the golden channels (i^e 



i^/x) |15l[ | was considered. It is found that the complementarity between these channels 
can be of great benefit if statistics are limited. There are also indications that the addition of these 
platinum channels is vital in order to resolve the degeneracy between standard oscillation parameters 
and non-standard interactions. Full details of the experiment can be found in Ref. 11[. In brief, 
we studied a setup ha ving a baseline of 1300 km with a beam capable of delivering 1.4 x 10^^ useful 
muon decays per year 1521 ] for 5 years per p olar ity. For the detector we considered either a 20 kton 
totally active scintillating detector (TASD) jl53l | and a 100 kton liquid argon (LAr) detector [154l |. 
both of which would be magnetised. These detectors would be capable of detecting and identifying 
the charges of both electrons and muons, providing access to multiple oscillation channels - the 
{v'^) disappearance channels, as well as the golden and platinum channels. 
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5.2.1. Sensitivity to standard oscillation parameters 



The 3<T CP discovery potential and sensitivity to the mass hierarchy of the LENF with either a LAr 
detector or a TASD is shown in Fig. [191 ft aken from Ref. 11 

1551], various /3-beams 



high energy neutrino fact ory ( HENF) 



). We also show t he sensit iviti es of the 



104 



113 



114 



1221 ]. T2HK [155|] and the 



wide-band beam (WBB) 156l ]. We see that a LENF with an optimistic LAr detector (left-hand edge 
of the blue band) has sensitivity to CP violation comparable to that of the HENF, for all values of ^13. 
A TASD also performs competitively for sin^(20i3) > 10~^. The sensitivity to the mass hierarchy is 
an order of magnitude better than that of the wide-band beam which uses the same 1300 km baseline. 
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Figure 19. 3cr a) CP discovery potential and b) hierarchy sensitivity of the LENF with a 20 kton 
TASD and 100 kton LAr detector, the HENF, /3-beanis, T2HK, and the wide-band beam (WBB). 



5.2.2. Sensitivity to non-standard interactions 



Non-standard interactions (NSI's) can be parameterised by the parameters Eq,^ [l57l . 11581 ]. which 
describe the rate of the transition I'a — ^ t^ fi- The g olden and platinum channels have leading order 
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1581 ]. We show the sensitiv ity o f the LENF with a 



1591 ]. simulated by MonteCUBES [l60|]. We illustrate how 



sensitivity to the parameters £e^ and Egr 
20 kton TASD to the NSI parameter 
the platinum channel enhances the sensitivity of the experiment, by showing the 68%, 90% and 95% 
allowed regions in the ^13 — plane, both with the platinum channel (red solid lines) and without 
(blue dashed li nes), for the case of = (Fig. [20b ) and Ee/j, = 0.01 (Fig. [20b). The current bound 



on is 0(1) [16l|. Fig. [20h shows that the LENF could improve upon this and Fig. 120b shows 



that Eefi = could be excluded at ~ 90% confidence if Ee^i ~ 0.01. We find that this sensitivity 
is not affected much by statistics, unlike for the case of standard oscillations. This is an indication 
that if NSI's are present, then the experimental sensitivity is limited by the degeneracies between 
the oscillation and NSI parameters, not by statistics. Unless this degeneracy can be resolved, for 
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instance by including the platinum channels, then the sensitivity to both standard oscillation and NSI 
parameters will be severely limited. 
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Figure 20. 68%, 90% and 95% allowed regions in the 6*13 



plane, with and without the platinum 



channel, for true values of a) = and b) 



0.01. 



5.3. The r-contamination of the golden muon sample at the Neutrino Factory jl2 | 



The problem of r-contamination was first studied in Ref. 1621 ] in the context of precision measurements 
of the atmosp heric parameters (Amgg, ^23) using the i^^ — )• channel at a Neutrino Factory. As it was 



shown in Ref. [I63l ]. it is useful to add all muons in the final state without muon charge identification, 

as the overall efficiency and threshold are much better and the backgrounds are still under control. The 

improved efficiency in the low-energy part of the neutrino spectrum, however, has the drawback that 

a previously irrelevant background becomes now potentially harmful. Oscillations into z^r's, otherwise 

1 



suppressed by low efficiency at low energy, enhance both the right- and wrong-sign muon sample 
Oscillations of z^e, — )■ will produce r's through i^t-N CC interactions within the detector that will, 
eventually, decay into muons (approximately 17% of them). These muons from taus will, therefore, 
"contaminate" the "direct" muon samples (coming from fg, — ?• i'^ oscillations). Notice that muons 
from taus are not background but as good signal as the direct muons. 

It is very hard to remove the muons from taus by applying kinematical cuts at an iron detector. 
Any cuts that attempt to do so drastically reduce the direct muon events as well and hence worsen 
the sensitivity to the oscillation parameters. They escape essentially all fi lters designed to kill the 
dominant backgrounds and directly add to the direct muon sample, see Ref. |l62| . On the other hand, 
neglect of the tau contribution will lead to an incorrect conclusion about the precision achievable at 
the Neutrino Factory on a given observable. The r -contamination must be added to the signal and it 
must be studied together with it. 

The problem of r -contamination of fg was studied in detail in Ref. [l^ . A good signal- 

to-background ratio is crucial to simultaneously determine with good precision ^13 and 6, since in 



They actually contaminate the electron sample, also 164 1 
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this channel we have a statistics of tens of events at most. To separate high-energy charged currents 
from the low-energy dominant neutral current background, is important a good reconstruction of 
the neutrino energy. For this reason, in the standard MIND analysis at the Neutrino Factory, the 
neutrino energy is reconstructed by adding the energy of the muon and that of the hadronic jet. This 
operation, however, yields a wrong result when the muon comes from a tau decay and it is detected 
at an iron calorimeter such as MIND, since no additional information regarding the neutrino missing 
energy in the r — t- Vr^p^^~ decay can be provideclfl. The sample of wrong-sign muons from the decay 
of wrong-sign taus will be distributed erroneously in neutrino energy bins, thus contaminating the 
wrong-sign muon sample by events whose parent neutrino energy is reconstructed wrongly. 

Consider a u-r of energy Ey^ , interacting in MIND and producing a wrongsign r of energy together 
with a hadronic jet of energy E^. After r -decay, Ey^ = E^ + E^ = {E^ + Emiss) + E^, where Emiss 
is the missing energy carried away by the two neutrinos in the r -decay. Experimentally, we observe 
the secondary muon and a hadronic jet, a signal essentially indistinguishable from that of a wrongsign 
muon from CC interactions. However, in this latter case, the addition of the (primary) muon energy 
E^j_ and of the hadronic jet energy E^ results in the correct parent energy, E^^ = E^ + E^. On the 
other hand, in the former case the addition of the (secondary) muon energy Ef^ and of the hadronic jet 
energy Eh results in the wrongly reconstructed fake neutrino energy Efake = + Eh = Ey^ + Emiss- 
If we divide the r three-body decay energy distribution in discrete fake neutrino energy bins, we find 
that for a monochromatic v^- beam of energy Ey^, the final muon will be assigned to a given fake 
neutrino energy bin of energy E^ with probability Vj{Ey^), where j = 1, ...,Nl^^. We can compute 
the distribution of z^r of a given energy Ey^ and divide them into I'r energy bins of energy EJ , where 
i = 1, N^in- '^^^ ensemble of the probability vectors Vj{El), for i and j running over all the and 
Vr energy bins, is represented by the migration matrix Mij , Fig. \n\ 

After having computed Mij , the number of total wrong-sign muons in a given neutrino energy bin 
is given by 



j=l,Nti„ 



(19) 



In Fig. [22]it is shown the fraction of muons coming from r-decay that can be found in the wro ng-s ign 
muon sampl^ after binning in the reconstructed neutrino energy, using the MIND efficiency 153] at 



L = 4000 Km (left) and L = 7500 Km (right). The data are shown for 6*13 = 2°, = ±90° (top and 
bottom, respectively) and the two possible neutrino mass hierarchies (red: normal hierarchy; blue: 
inverted hierarchy). In all cases, a significant r -contamination can be observed below 5 GeV, ranging 
from a minimum of 10% (inverted hierarchy, L = 4000 Km, 5 = 90°) to a maximum of 60% (normal 
hierarchy, L = 7500 Km, 6 = 90°). For normal hierarchy, the contamination in the energy range 
Ey S [5, 10] GeV drops to the percent level (with the only exception of L = 4000 Km, 6 = 90°). On 



^ This would not be the case at ECC [l65j or Liquid Argon detectors [l44l |. capable of separating the Ue signal 
from the Vf^ one and to measure precisely the kinematic of the process under study (so that the neutrino energy 

can be reconstructed precisely even — !> i^r — > t — !> ^ transitions). 

The setup is a 25 GeV Neutrino Factory with 5 x 10^" useful muon decays per year aiming at two 50 Kton MINDs 
located at L = 4000, 7500 Km, with 5 year of running time per muon polarity. An overall 2% systematic error has 
been considered. 
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Figure 21. The migration matrix Mij. Mij has been statis tical ly computed through a MonteCarlo 
simulation of the events, using the GENIE neutrino generator 166j with 10^ simulated z/rS per neutrino 
energy bin and 25 bins in the range E^^ £ [0, 25] GeV. From Ref. 12 1. 



the other hand, for inverted hierarchy the decrease of the r -contamination with the neutrino energy 
is softer, being as large as 30% for L = 7500 Km, S = 90°. Above 10 GeV, however, the contamination 
is at the percent level for both hierarchies (again, with the only exception of L = 4000 Km, 6 = 90° 
for inverted hierarchy). 

The T -contamination introduces, if not properly treated, an intolerable systematic error, in partic- 
ular for large 6*13. Fig. [23] (left), from Ref. [1^, shows that the test of the hypothesis that simulated 
data including the r -contamination, Ni{9i3,6), can be fitted using the direct wrong-sign muon dis- 
tribution, Nj'{9i3,5), fails at more than 3a for 6*13 > 5°. For 6*13 G [1°,5°], even if iVf (6*13,(5) can fit 
the r -contaminated data (with a relatively poor x^, though), the error in the joint measurement of 
^13 and 6 can be so large that it could actually prevent the use of the Neutrino Factory as a precision 
facility (see, again, Ref. [l3| for a detailed analysis of the errors introduced by a wrong treatment of 
the r -contamination). On the other hand, once Mij has been statistically computed, experimental 
data distributed in reconstructed neutrino energy bins, can be fitted using the complete wrong-sign 
muons distribution Ni{0i3,5), properly taking into account the r -contamination of the golden muon 
sample. Using this procedure, the systematic error introduced by the muons from taus is completely 
removed. The remaining error is the statistical error of the migration matrix elements, that is under 
reasonable control. 

It is worth noting that the r -contamination of the wrong-sign muon sample, once properly treated, 
does not worsen the measurement of ^13 and 6, as it was the case for the atmospheric parameters 
measurement. Fig. [^51f right) shows the comparison of the CP-fraction when only golden muons are 
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Figure 22. The fraction of r-contamination of the golden niuon sample as a function of the recon- 
structed neutrino energy, for = 2°. Left: L = 4000 Km; Right: L = 7500 Km. Top: S = 90°; 
Bottom: S =?90°. In red, we present the results obtained with normal hierarchy; in blue, with inverted 
hierarchy. 



considered (solid blue line) or when the total wrong-sign muon signal is taken into account (dashed red 
line) , as a function of sin^ 2^13 . It can be seen that the only difference between the two lines is a slight 
displacement of the wiggles at sin^ 2^13 10~'^(^?i3 ~ 1°)- The wiggles are a consequence of the loss of 
sensitivity to CP violation introduced by the so-called sign clones for negative 6 (a phenomenon known 
as vr-transit, 1671]). Since the location of the clones in the two samples differs, a small difference in the 



location of the wiggles is found when the two lines are compared. We can see, however, that once the 
r -contamination is properly treated, no (significant) loss in the CP-fraction is found everywhere else. 
This is a consequence of the extremely low statistical weight of the r -contamination in the golden 
channel. Since the signal is represented by tens of events, the r -contamination amounts to a few 
events, at most. Once the problem of the wrong assignment of muons from taus into reconstructed 
energy bins is solved by means of the migration matrix approach, the residual statistical impact of 
the r -contamination for ^13 < 10° is marginal. 
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Figure 23. Left: Test of the hypothesis that a simulation of the data that includes the effect of the 
T-contamination at L = 4 000 km can be fitted with the golden muon theoretical distribution. In the 
regions to the right of the contour lines, the hypothesis can be rejected at 1, 2 or Str (from left to right), 
assuming the goodness-of-fit statistics follows the distribution with n = 8 dof. From reference [1^ . 
Right: Comparison of the CP-fraction achievable at the IDS-NF baseline Neutrino Factory when the 
r-contamination is correctly taken into account (dashed red line) with respect to the ideal one in which 
no muons from taus are considered (solid blue line). 
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6. Physics beyond three-family oscillations 



The new facilities proposed within EUROnu have some potential to look for new physics beyond the 
standard three-family oscillations, such as Non-Standard Interactions, violation of unitarity of the 
PMNS matrix and light sterile neutrinos. In this section we present the WPG results on: neutrinoless 
double-beta decay 168]; sterile neutrinos as a solution to explain 2010 MiniBooNE data [3]; benefits 
of adding NSI's to sterile neutrinos to explain the 2010 MiniBooNE data [l^ : the potential of the 
Neutirno Factory to search for sterile neutrinos jisl ; the potential of the Neutrino Factory to distinguish 
NSI's from violation of unitarity of the PMNS matrix [16|]; and, NSI's in the Zee-Babu model 17 1. 



6.1. Neutrinoless double beta decay in seesaw models 



Neutrino oscillations, implying the massive nature of neutrinos, constitute an evidence for physics 
beyond the Standard Model (SM). Thus, models accommodating this neutrino masses become an 
important component in the search for new physics. In this context, one of the most promising 
processes is the neutrinoless double beta decay (0;//3/3 decay), whose detection is the purpose of several 
ongoing and upcoming experiments 70|. Since this proc ess is lepton number violating, its observation 
would imply that neutrinos are Majorana fermions 169l |. 

Among the i nost popu lar mode ls for neu trino masses, we find th e different types of seesaw mech- 

lese are the extensions of the 



183t | after the extra mediators 



anisms: type-I [170Ml73l ]. type-II |l74l4l78l ]. and type-Ill |l79l4l82| . T 
SM particle content that lead to the Weinberg d = 5 effective operator 
have been integrated out. 

The extra degrees of freedom associated to the just mentioned models, and required to induce the 
Majorana nature of the SM neutrinos, can also contribute to the Oz^/3/3 decay. The effects of the SM 
neutrinos and these extra states, introd uced to g enerate the neutrino masses, are usually analyzed as if 
they were independent (see, e.g., Refs. 



184 



188( 1). However, if the light neutrino masses are generated 



in the context of the model, important constraints between the light neutrino contribution and the 
one associated to the extra states arise. Here we will study the phenomenological implications of these 
correlations. 

In order to take into account the extra neutrino contributions to the Oz^/3/3 decay rate in the most 
general type-I seesaw scenarios, we have computed the nuclear matrix element (NME) without any 
assumption on the mass of the neutrinos mediating the process. They are related to the observable 
Oi^/3/3 decay rate as follows: 

2 

rOi//3/3 



In 2 



Goi 



(20) 



where Gqi is a well-known kinematical factor, nie the electron mass, mj the mass of the j neutrino, 
M^'^^^ (mj) the NME, which contain the dependence on the neutrino propagator, and Uej are elements 
of the neutrino mixing matrix. 
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The calculation was performed using Interacting Shell Model ( ISM ) nuclear wavefunctions, one of 
the most popular methods employed to obtain {)vj3j3 decay NMEs The results show two distinct 

regions for the behaviour of the NME as a function of the virtual neutrino mass: almost constant up to 
rrii ~ 100 MeV and then decreasing quadratically as the neutrino mass increases beyond 100 MeV. The 
transition mass 100 MeV is the typical momentum exchange of the decay, which is the momentum of 
the virtual neutrino |p|. This scale is fixed by the typical distance between the two decaying nucleons, 
r ~ 1 fm, and is sometimes referred as the nuclear scale^ 

The results of this computation are available at Ref. 190f | . Different approximations are required in 
order to perform the calculation, both obtaining the wavefunctions and in the treatment of the two- 
body transition operator, and consequently some uncertainties are induced into the NMEs. These can 
be estimated in % for light neutrino exchange {rrij < 100 MeV) and ~^4q ^ for heavy neutrinos 
{rrij > 100 MeV), more sensitive to the short range part of the transition operator. 

Within the type-I seesaw, depending on whether the extra mass eigenstates fall in the light or heavy 
neutrino mass regimes we can split their respective contributions to the amplitude: 

light light heavy 

Ao^Y, mU^iM^^^^irm) + ^ miU^jM°''^^{mi) + miU^jM'^''^^ (mi), (21) 

i I I 

where capital letters denote the mass index of the mostly sterile states and lowercase letters the mostly 
active SM states. Moreover, the diagonalization of the neutrino mass matrix provides an important 
constraint between the light and extra parameters: 

light light 

|]m,C/2, + |]m,C/2, = 0. (22) 

This simple relation, stemming from the fact that a Major ana mass coupling for the active neutrinos 
is forbidden by the gauge symmetry of the SM, provides a useful constraint that should be always 
satisfied. We can now distinguish three cases exhibiting very different phenomenologies depending on 
the mass regime of the extra mass eigenstates: 
AH extra mass states are light. In this case Eq. (|22p implies 

light 

A ^ - Y, miUli (mO'^^^(O) - MO^^^(m/)) . (23) 

Since in this regime the NMEs are basically independent of the neutrino mass, rrij = 

jVf 01^/3/3(0) ~ M^^/^f^imi) the rate of Oi^/3/3 decay is very suppressed. Indeed, only the different neu- 
trino masses in the NME prevent a full cancellation leading to a suppression driven by Am?/p^ with 
~ (100 MeV)2. Therefore, in this regime, the 0^(3(3 decay becomes experimentally inaccessible 
even being the neutrinos Majorana particles. 

All extra mass states in the heavy regime. Now the NMEs for this extra states are very 
suppressed compared to the SM ones. Furthermore, Eq. (j22p implies that 

heavy light 

A^-Y ^lUli (mO^^^(O) - M^^PP^mi)^ « Y rriiUlM'^^^^iQ). (24) 
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Figure 24. Bounds f rom CUORICINO on the extra neutrino mixing from QiyfS/S decay in ^^^Te, with 
a 90 % CL half-hfe 191|. We have assumed the extra neutrinos are non hierarchical and show the 
bounds as a function of their common mass. We compare the case in which the contribution from the 
SM neutrinos is properly taken into account (striped area) to that in which only the extra contribution 
is considered (above the red line). 



Thus, the contribution from the light active neutrinos (first term) dominates the transition rate. 
However, a much stronger bound on the mixing with the heavy neutrinos than the one usually shown 
in the literature can be extracted from it. This is shown in Fig. [24] using our calculation for the NME 
and the results from CUORICINO. 

Extra mass states in the light and heavy regimes. In this scenario the leading terms stem 
from the light states: 

light light 

^ oc ^ M°'^'5^(mi) + ^ miU^jM^"'^^ (mi). (25) 

i I 

Now it is possible to satisfy Eq. (j22p even in a situation where rriiU^- <^ mjU^j, by canceling the 
contribution of the extra heavy states against that of the extra light ones while keeping the light 
neutrino masses small. In such a situation, the contribution of the light extra states could dominate 
over that of the active. As an example, if we consider the controversial Heidelberg-Moscow claim 



for a positive decay signal [192], the accommodation of this signal through only SM neutrinos 

would require, using our ISM NMEs, 0.24 eV < ?ti^^ < 0.89 eV at 2a. The interpretation of this 
claim as light act ive SM n eutrinos is very disfavoured by the constraints from cosmology and neutrino 



oscillation data 



193l4l95l ]. However, this signal could be accommodated in a model with heavier 
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MB 


LS 


MB+LS 


KA 


MB+LS+KA 


Re 


(MB+LS+KA)+Re 


No Osc. 


A 


17.8 


15.0 




6.8 




51.0 






NDF 


16 


5 




8 




56 






GoF 


0.33 


0.010 




0.55 




0.66 




Osc. 


A mm 


10.6 


1.4 


13.3 


6.4 


24.1 


48.5 


75.6 




NDF 


14 


2 


18 


6 


26 


54 


82 




GoF 


0.71 


0.51 


0.77 


0.38 


0.57 


0.69 


0.68 




sin^ 2i9bf 


0.83 


0.0058 


0.0059 


0.0010 


0.89 


0.042 


0.013 






0.069 


8.13 


4.57 


6.76 


0.055 


1.86 


0.46 


PG 


AXmin 






1.30 




5.80 




2.94 




NDF 






2 




4 




2 




GoF 






0.52 




0.21 




0.23 



Table X. Values of x^i number of degrees of freedom (NDF), goodness-of-fit (GoF) and best-fit 
values sin^ 2i?bf, ^"T-m of the oscillation parameters obtained from the fit of various combinations of 
MiniBooNE (MB), LSND (LS), KARMEN (KA) and reactor Bugey and Chooz (Re) antineutrino 
data. The first three lines correspond to the case of no oscillations (No Osc). The following five lines 
correspond to th e case — >■ De oscillations (Osc). The last three lines give the parameter goodness- 
of-fit (PG) [l97|. The variations of sin^ 2i9bf and Am^f depending on the fitted data sets are due to 
the oscillating character of Pp^_>p^ in Eq. (^5)) . See Ref. 



131 for details. 



neutrinos (which are not bounded by cosmology) mediating the process. Indeed, we could reinterpret 
the result as 0.24 eV < ^^cavy ^^^2^ ^ ^ gg 

As for the type-II and type-Ill seesaws, cur rent bounds from accelerator experiments place the 
extra degrees of freedom in the heavy regime [196l |. This effectively reduces the situation to that 
which appears for the type-I seesaw with only heavy extra states. In the same manner, mixed type-I 
and type-II/type-III seesaw models (with extra neutrinos below 100 MeV) resemble the situation of 
type-I seesaw with both light and heavy extra states. 



6.2. Short-Baseline i/^ — > z/g Oscillations 13 | 



The MiniBooNE collaboration jl98l | recently reported the observation of a si. 



gnal of short-baseline 



z?^ —7- Ue transitions compatible with that observed in the LSND experiment [199l |. The agreement 
of the MiniBooNE and LSND signals in favor of neutrino oscillations is remarkable, because the two 
experiments observed the signal of — )• z^e transitions at different source-detector distances and 
different neutrino energy ranges. Since only the ratio of distance and energy is sim ilar i n th e two 



experiments and neutrino oscillations depend just on this ratio (see Refs. |4ll. l55l. ISg. |200I4204| ]). the 
neutrino oscillation explanation of the two signals is strongly favored. On the other han d, the Mini- 
BooNE collaboration did not observe any signal of short-baseline — ?• Ue transitions j205| compatible 
with the MiniBooNE and LSND signals of z/^ — t- z/g transitions. Therefore, it is possible that the 
effective parameters which govern neutrino and antineutrino oscillations are different, maybe because 
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of a violation of the CPT symmetry 206l4224i | . From a phenomenological point of view it is inter- 
esting to consider the neutrino and antineutrino sectors i ndepend ently, especially in view of possible 
experimental checks of the short-baseline — )• Ve signal [225l4228t |. Here we adopt this point of view 
13t | and we present the results of a combined fit of the MiniBooNE and LSND antineutrino data 
in favor of short-baseline V n — ?• v^. transitions, together with the constraints imposed by the data of 
the KARMEN experiment 229l | in which the transitions have not been observed. We also take into 



account the constraints imposed by the absence of short-baseline Ve disappearance observed in reactor 
antineutrino experiments. 

In the analysis of the data of z/^ — )■ oscillation experiments we consider the simplest case of an 
effective two-neutrino-like short-baselin e oscillati on probability, similar to that obtained in the case of 
four-neutrino mixing (see Refs. 
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5S 



20C 



2031]) 



P. ^p^{L/E) =sm^ 2^9sin2 



Am^L 
AE 



(26) 



where Am^ is the relevant neutrino squared-mass difference and •& is the effective mixing angle for 



— )■ transitions. 



Tab. |X] shows the results of the fit of various combinations of MiniBooNE (MB), LSND (LS), 
KARMEN (KA) and reactor Bugey 230 ] and Chooz [s^] (Re) antineutrino data. 
The left panel in Fig. [25] shows a superposition of the 90% and 99% C.L. allowed regions in the 



sin^ 2'&-AmP' plane obtained from the combined fit of MiniBooNE, LSND and KARMEN Un — )• z/g data 



and the exclusion curves obtained in Ref . |231 1 from the fit of reactor Bugey and Chooz z/g — ?• I'e data, 
which currently provide the most stringent constraints on short-baseline reactor z^g disappearance. 
The mo del- independent inequality 



(27) 



implies that in the left panel in Fig.[25]the large-sin^ 2-!? part of the straight region below Am? ~ 2eV^ 
allowed by the combined fit of MiniBooNE, LSND and KARMEN z^^ — ^ z/g data is excluded by the 
results of reactor antineutrino experiments. Quantitatively, only the parts with sin^ 2'& < 3 x 10~^ 
and sin^ 2'!9 < 5 x 10~^ are allowed at 90% and 99% C.L., respectively. 

The right panel in Fig. [25] and the last column of Tab. [X] give the results of the combined fit of 
accelerator MiniBooNE, LSND and KARMEN i)^ — )■ z/g data and reactor Bugey and Chooz — )• z/g 



data assuming an equality in Eq. (j27p . The value of the parameter goodness-of-fit in Tab. [XI shows 
that the accelerator and reactor data are compatible under the hypothesis of — )• z/g oscillations. 

From the right panel in Fig. [25] one can see that there is a favorite region at about 95% C.L. around 
the best-fit point for 2 x 10~^ < sin^ 219 < 5 x 10~^ and 0.2 < Am^ < 2eV^. Larger values of Am^ 
are allowed only at more than about 95% C.L. for 5 x 10^^ < sin^ 2-;? < 5 x 10"'^. 

This region is interesting for a study of the possibilities to check t he LSND and MiniBooNE indica- 



tion of short-baseline — t- z^g oscillations with future experiments 225142281 1. 
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Figure 25. Left: superposition of the allowed regions in the sin^ 2'd-Am^ plane obtained from the 
combined fit of MiniBooNE (MB), LSND (LS) and KARMEN (KA) De data and the exclusion 

curves obtained from the fit of reactor Bugey (Bu) and Chooz (Ch) data. Right: allowed 

regions in the sin^ 2-d-/S.m?' plane and marginal A^^ — ~ Xmin'^ fo'^ ^i^^ 2i9 and /S.rn? obtained 
from the combined fit of MiniBooNE (MB), LSND (LS) and KARMEN (KA) D^, data and the 

exclusion curves obtained from the fit of reactor Bugey and Chooz (Re) Vf, — > Vf. data. The best-fit 
point is indicated by a cross. See Ref. for details. 



6.3. MiniBooNE/LSND data: NSI's in a (3+l)-scheme vs. (3+2)-oscmations []J] 

Recent MiniBooNE anti-neutrino data 



LSND evidence 



1981 ] indicate an excess of Ve events, in agreement with the 



199( 1 for — )• z/e transitions. It is known that oscillations with one or more sterile 



neutrino(s) at the eV scale are not sufficient to explain the global data 232, 12331 ] . The preprint 



EURONU-WP6-10-24, Ref. Q, investigated the possibility that in addition to a sterile neutrino there 
are some non-standard neutrino interactions (NSI), beyond the Standard Model weak interactions. 
Since matter effects are tiny for the short baselines relevant here, charged-current (CC) type NSI in 
the neutrino source and detector are considered. Thanks to the interference between NSI effects and 
oscillations with Am^^ ~ 1 eV^ CP violation is obtained, even in the presence of only one mass scale. 
This effect is used to reconcile the indication for i/^ — )• Uf. in anti-neutrino experiments (LSND and 
MiniBooNE) with the absence of a signal in MiniBooNE neutrino data. 

A general parameterisation of the relevant transition and survival probabilities in the presence of 
oscillations (within the one- mass scale approximation) and NSI is presented, and particular com- 
binations of mixing matrix elements UaA and NSI parameters £ap entering in the probabilities are 
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identified. This drastically reduces the number of independent parameters and makes a general fit to 
global short-baseline data feasible. 

Two versions of the (3+1) NSI model are considered. In the general case (denoted NSI^) one 
makes use of the fact that the neutrino production mechanism in LSND (and in KARMEN) is muon 
decay (purely leptonic), whereas in all other experiments neutrino production and detection are semi- 
leptonic, involving transitions between u and d quarks. Therefore, in the presence of suitable NSI 
parameters one can decouple the transition probabilities in LSND and KARMEN from the rest of the 
data. In this case one obtains an excellent fit to the global data and the tension between appearance 
and disappearance experiments is resolved. In this case MiniBooNE does not provide a direct test of 
LSND, since different combinations of parameters are relevant for them. Also, in the global fit the 
excess observed in MiniBooNE anti-neutrino data is not reproduced. 

For the second version of the (3+1) NSI model the assumption is adopted that NSI involving the 
charged muon can be neglected. In this case exactly the same NSI parameters are relevant for LSND 
and KARMEN as for all other experiments. In this constrained model (NSF) one makes use of the CP 
violation due to NSI-oscillation interference to reconcile neutrino and anti-neutrino data. It is shown 
that in the NSF model there is a factorisation between appearance and disappearance amplitudes, 
similar to that in the (3+1) oscillation scheme. Therefore, it is more difficult to satisfy constraints 
from disappearance experiments and some tension is left in the fit. However, also this model provides 
significant improvement of the global fit compared to the pure oscillation case. 

The results of the fits are presented in terms of effective parameters, representing the specific combi- 
nations of NSI parameters entering in the transition probabilities. However, for both cases, NSF and 
NSP, examples are provided of how to realise the required parameters in terms of the fundamental 
mixing and NSI parameters. It is shown that values in safe agreement with bounds on the various e's 
can be found to realise the fits which require e's of order a few xlO~^. 

The quality of the (3+1) NSI fits are compared to an updated fit in the (3+2) oscillation scheme, 
which also allows for CP violation due to the presence of two relevant mass scales. Similarly to (3 +1) 
NSI, in (3+2) the appearance experiments can be described very well. However, as previously 2321 ] . for 
(3+2) oscillations significant tension remains in the global fit between appearance and disappearance 
experiments. The improvement of (3+2) compared to (3+1) is not significant, in terms of gain 
per new parameter. None of the scenarios considered in [ij] can explain the MiniBooNE low energy 
excess of events when disappearance data are taken into account. Therefore, the data below 475 MeV 
is excluded from the analysis, relying on a separate explanation for this anomaly. 

The predictions of the model for future experiments depend on the detailed realization in terms 
of mixing and NSI parameters. In general one may expect some signals in searches for deviations 
from the standard three-fiavour oscillation picture in both respects, sterile neutrino oscillations as 
well as NSI. Several proposal s to sear c h fo r ste r ile n eutri nos a t the eV scale have been presented 



recently, see for example 15|, l22ll . |225| . 12271 . |234| . |235| | . In 236l | implications of sterile neutrinos for 



latest cosmological data have been investigated. Recen t studies on N SI in the context of upcoming 



and far future experiments can be found, e.g., in [16|, Il57l . 12371 . l238l |. A specific prediction of the 
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scenario are zero-distance effects in appearance searches jl6l. l239l . |240| ]. Hence, the observation of an 
energy independent appearance probabihty at very short distances is a characteristic signature from 
this kind of models. 

The model may also provide a signature at the LHC. Typically, realising CC-like interactions as 
the ones considered here require a charged particle as mediator. The NSI parameters e measure 
the strength of the new interactions relative to the standard weak interaction strength set by Gi;'. 
Therefore, from the fit results, e ~ 0.01, one expects that the mass of a mediator for a dimension- 
6 operator should be roughly one order of magnitude larger than the W boson mass. Hence, one 
might expect charged particles to s how up at the TeV scale, with good prospects to be observed at 
LHC. However, the results of 



241 



242l | suggest that NSI at the level of 0.01 are difficult to obtain 



from dimension-6 operators without being in conflict with bounds on charged-lepton processes. As 
discussed there, a possibility to obtain such large NSI would be to go to dimension-8 operators and 
allow for some fine tuning. 



6.4. Sterile neutrinos beyond LSND at the Neutrino Factory 15 



In this section we will discuss the discovery reach for sterile neutrinos at the standard IDS-NF baseline. 
We consider the simplest scenario in which only one sterile neutrino is added to the three active 
neutrinos of the Standard Model in the so- called 3+1 sc heme which recovers the standard picture in 
the case of small active-sterile mixings [15l . 



234 



243l4248t | . For the sake of simplicity, we just focus on 



the scheme where the fourth state is the heaviest and normal hierarchy is assumed in the standard 
sector. The numerical results on the discovery reach can be understood from the analytical expressions 
of the transition probabilities; we found particularly illuminating the following parameterization of 
the 4x4 unitary mixing matrix U: 



U — i?34(^34, 0) i?24(^24, 0) i?i4(^i4, 0) i?23(^23i ^s) ^13(^13) ^2) -^12(^12, ^1 



(28) 



where Rij{9ij, 61) are the complex rotation matrices in the ij-plane. In the short-baseline limit 
IA41I = Am1iL/AE ~ 0{1) ^ |^3i|) the matter effects can be safely ignored, and the relevant 
probabilities read: 



Ve^i = Vfie = 4:cl^sl^sl4^ sin^ A41 
Pee = 1 - sin^ (26ii4) sin^ A41 

'P^.T = 4Ci4C24S24«34 sin^ A41 

Vf,f, = l- c\^s\^ [3 + 2cf4 COS (2024) - COS (2^14)] sin^ A41 



(29) 
(30) 
(31) 
(32) 



where we used the short-hand notation Ajj = Am? L/4£'. From these probabilities we can see that 



024 can be measured by V^^ and ^14 by Vee- On the other hand, ^34 only shows up in combination 
with the other small mixing angles. For long baselines, some of the relevant features of the probability 
transitions can be well understood using simple perturbative expansions: for A31 = 0(1) ^ A41 and 
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Figure 26. The exclusion limit at 90% CL (2 d.o.f) for sin^ 2&i4-ATO|j^(i ~ 1,2,3) (region on r.h.s. 
of curves excluded). Here gives results of the standard IDS-NF (AOOO km and 7500 km) setup with 
near detectors ( thick solid curves ), and the impact of the near detectors is also shown separately ( thick 
dashed curves). The thin dashed curves only consider systematics. 



up to the second order in S13, S14, S24, 534, S23 = sin ^23 — and considering A21 as small as s?^-, we 
obtain: 



V^,^, = cos^(A3i)(l - 2s^4) + 8s^3 sin^(A3i) + €{2^12 sin(2A3i) + 

2S24S34 cos (53 A„ sin(2A3i) - (33) 

.2a ^ a ^ (A31-Ae)AeSin(A3i) - A31 sin (A31 - Ae) sin(Ae) 
^■Sl3^31 COS(/i3ij 



(A31 - Ae)2 

V^,r = sin2(A3i)(l - 85^3 - - s^j _ ^^/^^^ sin(2A3i) - 

S24S34sin(2A3i) [2A„cos(53 - sin53] - (34) 
A31 {sin(A3i - Ae) + sin(Ae)} - 2(A3i - Ae) Ae cos(A3i) 



S13A31 sin A31 



(A31 - A 



from which we learn that at the long baselines, ^24 is best accessible by V^^ with the first term 
proportional to cos^( Ai3). The leading sensitivity to ^34 can be expected from Vfj^r (the discovery 



channel as claimed by [23j]). Notice also the dependence on the phase 5-^ in both probabilities which 
makes them useful to check whether other sources of CP violation beside the standard one (^2 in our 
parameterization) can be tested at the Neutrino Factory. Let us first discuss general constraints to 
the new mixing angles ^m, ^24) and 6*34 and the additional mass squared difference /^m\^ without any 



additional assumptions [15[. The considered oscillation channels, simulated with a properly modified 
version of the GLoBES software, are electron to muon neutrino (appearance channels) and muon 
to muon neutrino (disappearance channels) oscillations. Beside the standard IDS-NF setup, we also 
considered near detectors with a fiducial mass of 32 1 and a distance of d = 2 km from the end of the 
decay straight, which corresponds to the effective baseline of 2.28 km. We show three different exclusion 
planes Oi^-Am^i in Fig l6.4i The main sensitivity is obtained at about Am|^ ~ 10 eV^, which comes 
from the distance chosen for the near detectors. Since the efficiencies for muon neutrino detection 
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are typically better, the sensitivity to 624 is slightly better than that to 614 for our assumptions. As 
expected, there is no sensitivity to 634 coming from the near detectors, because the i^r disappearance 
channel does not exist. For the effect at the long baselines, it is first of all useful to consider the thin 
dashed curves with systematics only. In all three panels, the sensitivity changes as a function of Amj^^ 
in the region where Amli ~ Am'^i. It comes from the fact that the Neutrino Factory is sensitive 
to the atmospheric oscillation frequency, whereas for Amli ~ Arn^i, no particular additional effects 
from the solar frequency can be found. As expected (see Eq. ()33p ). the main sensitivity is found for 
624. However, there is also some sensitivity to 9i4, which vanishes after the marginalization, and some 
sensitivity to ^34, which is even present for Am|^ = for systematics only. After marginalization (thick 
solid curves), only the sensitivities to 624 and ^34 remain in the Am|^ regions close to the atmospheric 
Am|^ and above, where the effects of Am|-^ average out. Very interestingly, note that mixing angle 
correlations destroy the sensitivities for Am^i = Arri^i, where m4 = and no additional Amli 
is observable, leading to small gaps (see horizontal lines). We have tested that the sensitivity to 
^34 is a matter potential-driven, statistic limited higher order effects present in the muon neutrino 
disappearance channels. In view of the three panels, it is not easy to disentangle the parameters for 
arbitrary massive sterile neutrinos. Parameter correlations lead to a pollution of the exclusion limit 
of a particular mixing angle with Am|^. In addition, there is a competition between Amli and Am^i 
at the long baseline. Near detectors, on the other hand, have very good sensitivities to 614 and 624 
but cannot measure ^34. Nevertheless, the absolute values of the sensitivities are quite impressive. 
Especially, 624 can be very well constrained close to the atmospheric mass squared difference range. 
This indicates that sterile neutrino bounds in that range should be also obtainable from current 
atmospheric neutrino oscillation experiments. 



6.5. NSI's vs. non-unitary lepton flavor mixing at a neutrino factory 16 | 



Apart from the study of neutrino oscillation parameters, one of the most important tasks in future 
neutrino oscillation experiments, e.g., a neutrino factory (NF), is to search for non-standard effects 
stemming from the exchange of heavy fields presented in high-energy theories. In the language of 
the effective theory, the impacts of heavy particles can be parametrized in terms of a tower of higher 
dimensional non-renormalizable operators, e.g., the dimension-five operator O^^ = (L^^ir^^) {(pir'^L) 
responsible for the generation of neutrino masses. According to the nature of the heavy fields, the 
dimension-six operators can be classified into two categories, the scalar me diated op erators and the 
fermion mediated operators (See discussions and references in, e.g., Refs. [itI. |249 . 250] on the possible 



theories behind these operators). These dimension-six operators typically suffer from severe exper- 
imental constraints on the lepton flavor violating processes, since i — >• Hi decays may be induced 
at tree level. However, there are two leptonic dimension-six operators which do not lead to charged 
lepton flavor violation, 

= {17- L) {L ■ L') , = (lJ) i^ U^L^ , (35) 
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where L and (j) stand for the Standard Model lepton doublets and the Higgs field. The scalar medi- 
ated operator O'^ leads to non-standard neutrino interactions (NSIs), whereas the fermion mediated 
operator O-^ results in non-unitarity (NU) effects. Note that, the NU effects are similar to these from 
NSIs, and therefore, both NSIs and NU can be described by using the same parametrization, i.e., the 
conventional parameters, where A = s,m,d correspond to non-standard effects in the neutrino 
production, propagation, and detection processes, respectively. 

The generic bounds on above two dimension-six operators come from rare lepton decays ia — ^/37, 
the invisible de cay width of the Z-boson and the universality test of weak interactions, and can be 



found in Refs. 
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24l[. An important question to be answered is if we can distinguish these two 



operators at future neutrino oscillation experiments so as to find hints on the origin of non-standard 
effects. To this end, we remark that the NU effects are fundamental and process-independent, whereas 
the NSI effects are experiment-dependent. For example, O'^ does not affect the neutrino production 
process at a superbeam experiment since it is essentially a hadron decay process. We summarize in 
Tab. IXII which non-standard effects are allowed for a neutrino factory and a superbeam if the origin 
are the discussed leptonic dimension-six operators. In addition, there exist non-trivial correlations 
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Tabic XI. Allowed parameters 
from the discussed dimension- 
six effective operator classes in a 
neutrino factory (v-factory) and 
a superbeam experiment (SB). 



between the source and matter non-standard effects for both O'^ and O-^ 242], i.e., 
ie^Jr (NSIs), e;r. = -e^. (NU), 



'tJ.T 



(36) 



which motivate us to investigate the possibility of distinguishing the origin of non-standard effects at 
a neutrino factory alone [16 1. 

In the left column of Fig. [271 we illustrate the potential of discriminating 0£ from O-^ at an IDS- 
NF (International Design Study) with parent muon energy = 25 GeV [2511 ] (See also a detailed 



description of the experiment setups in Ref. [3l). In our numerical analysis, we make use of the 
"true" parameters in the case of {O^), and then fit the data with only (O^). For the IDS-NF 
neutrino factory combined with several different near detectors, the curves in the left column show 
that there is just a very small region beyond the bound at the 90 % C.L. on the O"^ operators, 
where the data generated with O'^ can be distinguished from the O'^ even if the OPERA-like near 
detector at the longer baseline is used. If O'^ is simulated, however, it can be distinguished from O-^ 
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True value of |e' | True value of |e' | True value of |e'J 




ND-L: Large size (OPERA-like) 

fiducial mass 2 kt, d = 1 km 

ND-M; Medium size (e.g., SciBar-sized) 
fiducial mass 25 t, d = 80 m 

ND-S: Small size (e.g., silicon vertex-sized) 
fiducial mass 100 kg, d = 80 m 

OND@130km: Large size (OPERA-like) 
fiducial mass 2 kt, L = 130 km 



True value of \ True value of |f:^ | 



Figure 27. Left column: regions in the (]e^^|-(/)^^J-pZane where the simulated e^^ induced by one 
type of operator can be uniquely established, i.e., the other type of operator is excluded at the 90 % 
C.L. (regions on the right-hand side of the curves) at an IDS-NF (E^^ = 25 GeVj. The discovery 
reach is also displayed in each plot. Middle column: 90 % C.L. of CP discovery potentials in the NU 
framework (upper plot) and N SI framework (lower panel) at an IDS-NF. Right column: same plot as 
the lower panel of the left column, but for a low-energy neutrino factory (E^ ~ 4.5 GcV) alone and in 
combination with the ND-S and ND-S-hND-L. Taken from fl^ l- 



for a part of the parameter space beyond the current bound with ND-L. In the middle column of 
Fig. [28| we plot the CP discovery potential for both O'^ (upper panel) and O'^ (lower plot) induced 
CP violations. This is defined as the ensemble of true values of </)^^ (phase of e^,-), which cannot be 
fitted with the CP-conserving values = 0, zbvr at 90 % C.L. The combination of the standard IDS- 
NF neutrino factory with different large enough near detectors may discover CP violation, somewhat 
beyond the current bounds, especially for (^*^ ~ ±tt/2. There are no qualitative differences between 
and O^. Finally, we sh ow, in the right upper panel of Fig. [28| that a low-energy neutrino factory 
{E^ = 4.5 GeV) [l47|, [l49| is not able to exclude the operator at 90 % C.L., since the exclusion 
regions are excluded by current limits already. 

In conclusion, differentiating between NSIs and NU should be one of the key priorities of searches 
for new physics effects, since the nature of non-standard effect points towards the nature of the heavy 
mediator. The components necessary for this search are Vr detection at least in near detectors, both 
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at high-intensity superbeams and a neutrino factory. For the neutrino factory, a high enough muon 
energy is mandatory for the discussed non-standard effects searches, which means that the high-energy 
neutrino factory should at least be an upgrade option even for large ^13. In addition, for non-standard 
effect searches, the size of the near detector is very important, which means that for all applications, 
large enough detectors are needed. 



6.6. Non-standard neutrino interactions in the Zee— Babu model 17 



While the origin of neutrino masses emerges as one of the unsolved and fundamental problems in 
particle physics, the mechanism of radiative mass generation provides a natural method to obtain 
small neutrino masses. In such a framework, neutrino masses are exactly vanishing at tree level, 
and are in duced as finite radiative corrections. The simplest version of this kind of models, i.e., the 
Zee model [252l |. cannot accommodate current experimental da ta, since the predicted leptonic mixing 



angle 612 is too large. Alternatively, in the Zee-Babu model j253H255| |. two SU(2)i^ singlet scalars 
(one singly charged and one doubly charged h~^+) are introduced besides the Standard Model 
(SM) particle content, and neutrino masses can be generated at two-loop level. In addition to the 
generation of neutrino masses, the exchange of heavy scalars also results in lepton flavor violating 
processes such as /x — t- 3e and fi — )• 67, which can be dramatically enhanced compared to those in the 
SM. Most interestingly, if kinematically accessible, the new scalars could be directly produced at the 
Large Hadron Collider or the future International Linear Collider, and the decay of the doubly charged 
Higgs may induce very clean like-sign bi-lepton events, which are basically free from SM background, 
and opening the possibility of collider tests of the neutrino mass generation mechanism. 

Besides colliders, the next generation neutrino oscillation experiments, e.g., neutrino factories, will 
also help us to unveil the underlying physics behind neutrino masses, in particular, throug h exp loring 
non-standard neutrino interactions (NSIs). Unfortunately, in most of the low-scale models [249|], NSIs 
are typically lin ked t o lepton flavor violating processes of charged leptons, yielding too tight bounds, 



see, e.g., Refs. |161 



241 



2561 ]. However, in the case of the Zee-Babu model, the situation is more 



involved, since the masses of singly and doubly charged Higgs in principle can be well separated and 
a different set of Yukawa couplings controls charged lepton and neutrino interactions with the scalars. 
In the following we summarize the results of Ref. 
Zee-Babu model has been performed. 



171], where a detailed investigation of NSIs in the 



The Lagrangian of the minimal Zee-Babu model reads 



= jCsm + fai3LiaCi(T2LLi3h^ + ga/se^e/sk'^'^ - fih h k'^^ + h.c. + Vh , 



(37) 



where Ll denote left-handed lepton doublets, e are the right-handed charged leptons, and the scalar 
potential Vh contains the couplings among scalar fields. / and g are antisymmetric and symmetric 
Yukawa couplings, respectively. Neutrino masses are generated at two-loop level as 

1 



487r2 M2 



ifD,g^D,f, 



(38) 
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where D^, = diag(me, m^, m,-) contains thecharged-lepton masses, M = max(mfc, ?7t,/j), and I{r) is 
a dimensionless function of order unity j257| |. Furthermore, the tree-level exchange of induces 
non-standard lepton interactions via dimension-six operators 

£NSi ^ [l^PLVp) {v^PLe;) = 2V2GFe':, {V^rPLyp) [Tp^.PlL) , (39) 



mi 



where e^*^ ~ ^-^^ fa/sfpa { t^) ^ canonical NSI parameters. Using the conventional notation 



in the literature, e.g., Ref. j238l |. we find that the following NSI parameters are induced 



.m ee _ /e/3 fea s e/i _ f^iefer s e^i _ f IJ-tI efi 



where the superscript m and s correspond to NSI effects in matter and in the neutrino source at a 
Neutrino Factory, respectively. 

Equations ()38p and (|40p indicate that there exists a firm connection between neutrino parameters 
and NSI parameters. For example, in the of normal mass hierarchy (mi <^ m2 ^ m^, NH), the relation 
fe^i — fer — furf^ approximately holds, while in the case of inverted mass hierarchy (m2 > mi ^ ms, 
IH), l/e^l ~ |/er| and |/^T-| ~ |/eT|'Si3/s23 Can be obtai ned. Taking into account the experimental 
constraints from the universality in £a — )• itT^i^ decays [2571 ] . we can roughly estimate that \fe^\ ~ 
|/er| ^ 0.05 (m/j/TeV) and \f^r\ ^ 0.1 (m/i/TeV). Therefore, compared with Eq. ([1CT|) . the important 
NSI parameters in the NH case is e^^, while in the IH e!^, and e^^ can be sizable. 

Taking into account the experimental bounds summarized in Ref. [257l |. we present, in Fig. 1281 the 
allowed regions of NSI parameters and sin ^13 at Icr, 2a, and 3cj C.L. We find, from the left column, 
that and can only be sizable in the IH case, while the present long-baseline experiments are 
not very sensitive to these parameters. As for the source related NSI parameters, may reach 
values up to 10""^ for both hierarchies, whereas e^^ can be as large as few xlO~^ only in the IH case. 
In particular, for a scalar mass scale of 1 TeV, a non-trivial lower bound on the NSI parameters of 
order 10~^ is found, indicating that the model is rather constrained from the requirement of a correct 
neutrino mass matrix. The right column in Fig. [28] shows correlations between NSI parameters and 
the mixing angle ^13. Especially, in the IH case, one obtains a quite strong prediction for ^13, i.e., 
a lower bound sin^ ^13 > 10~^. Such a sizable lower bound is of particular interest, since it would 
guarantee a discovery at the forthcoming reactor or long-baseline experiments in the near future [l5Cll | . 

In conclusion, sizable NSIs can be accommodated in the Zee-Babu model in a natural way, and 
they depend on both the mass scale of the charged scalars and on the type of neutrino mass hierarchy, 
NH or IH. We further remark that, in addition to the NSIs, the singly and doubly charged scalars 
could be directly produced through the s-channel processes at the Tevatron and the LHC. Therefore, 
the interplay of the phenomenology at colliders, search for lepton flavor violation, and NSI effects at 
neutrino oscillation experiments could play a complementary role towards the goal of identifying the 
true mechanism of neutrino mass generation. 
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Figure 28. Allowed regions of NSI parameters and sin6'i3 at Icr, 2a, and 3(t C.L. We take m/,. = mfe = 
fj, = 10 TeV for the figures in the left column and irih = = fi = 1 TcV for the figures in the middle 
and right column. Taken from [17 1. 
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7. Physics potential of EUROnu facilities as of April 2011 



18| 



In Figs. [29ll3T] we present the comparison of the physics potential of the three EUROnu facihties 
(the 25 GeV Neutrino Factory, the 7 = 100 Beta-Beam and the SPL Super-Beam) as shown at the 
Technical Review Meeting in Geneva, April 13th, 2011. Results for the 7 = 100 beta-beam are shown 
for several values of the neutrino fluxes, as well as in combination with the SPL. In addition to the 25 
GeV Neutrino Factory, results for the 10 GeV Low-Energy Neutrino Factory are also shown in case 
sin^ 2013 were found to be larger than 10"^ (as it is indeed the case according to the recent results of 
T2K, MINOS, Daya Bay, RENO and Double Chooz experiments). 

The three figures correspond to: (1) sensitivity to ^13, Fig. [29l (2) CP discovery potential. Fig. [30l 
and, (3) sensitivity to the neutrino mass hierarchy, Fig. 1311 The curves correspond to 3a CL (1 dof). 
The other parameters have been fixed to ^m'^i = 0.0024eV^, Amg^ = 8 x 10~^ eV^ and 623 = 45°. 
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Figurc 29. Comparison of the sensitivity to 6*13 of different future facilities as a function of sin 2613. 
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Figure 30. Comparison of the CP discovery potential of different future facilities as a function 
of sin^20i3. Prep ared by P . Huber for the EUROnu Tech nical Rev iew, Apri l 13t h, 20 11, usin g th e 
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Figure 31. Comparison of the sensitivity to the neutrino mass hierarchy of different future facilities 
as a function of sin^ 26*13. Pr epared by P. Huber for the EUROnu Techn ical Revie w, A pril 1 3th, 2011 , 
using the GLoBES package [253, [253. Curves are taken from [a] ^0$, [h] , [c] [123, [d] [150 1 
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8. Summary of the NuFlavour workshop jl9( | 



In 8-10 June 2009, the workshop "Flavour physics in the era of precision neutrino experiments" at 
Coseners House, Abingdon, UK, focussed on a critical review of the physics case for neutrino physics 
and long baseline neutrino oscillations from a theoretical perspective. This workshop was organised 
with partial support from the European Community under the European Commission Framework 
Programme 7 Design Studies: EUROnu (Project Number 212372). The topics which were discussed 
included: i) lepton flavor violation (LFV) from grand unified theory (GUT) see-saw models and 
from TeV see-saw models, ii) Leptogenesis in the context of neutrino mass models, iii) A theoretical 
perspective on lepton flavor physics at the TeV scale within supersymmetric (SUSY) models and extra- 
dimension models, iv) Interplay between neutrino masses and other phenomenological signatures, v) 
Neutrino physics and cosmology /astroparticle physics, and finally vi) Performance indicators in long 
baseline neutrino oscillation experiments. 

8.1. LFV from GUT see-saw models and from TeV see-saw models 

See-saw models provide a natural solution to the smallness of neutrino masses and they can be em- 
bedded in Grand Unified Theories (talk by G. Ross) or emerge at the TeV scale, being testable at 
the LHC in the forecoming future (talk by E. Ma). All these models need also to account for the 
low energy flavor structure (talk by M. Schmidt). A series of theoretical questions emerged at the 
meeting: Is it possible to discriminate between GUT and TeV-scale see-saw models? How do these 
models inscribe into a more general theory? Where does the smallness of neutrino masses come from, 
and does it require fine-tuning anywhere? How can flavor mixing be implemented? Is it possible 
to define some "benchmark" scenarios? If so, which ones? Is it possible to have arguments for the 
precisions required for ^13, 5, and the Majorana phases? What information on the origin of neutrino 
masses can be obtained from the mass hierarchy measurement? These questions were addressed in the 
discussion session lead by M. Frigerio. It was mainly emphasized that there are two aspects which need 
to be studied: the origin of neutrino masses and the rationale for the leptonic fiavor structure. Three 
contributed talks were also given by A. Adulpravitchai on Non-abelian Discrete Flavor Symmetries 
from T'^ /Zj\f Orbifolds, by D. Hernandez on the Minimal Flavor Seesaw Models and by M. Satriawan 
on Possible mechanism for generating a very small Dirac neutrino mass. 

8.2. Neutrino physics and the cosmology /astroparticle physics complementarity 

The bath of relic cosmological neutrinos, which decoupled at T 1 MeV, affects the evolution of 
the Universe, in particular Big Bang Nucleosynthesis, large scale structure formation and the Cosmic 
Microwave Background and allows to constrain neutrino masses, the number of neutrinos and the 
type of neutrino interactions (talk by S. Hannestad). Sterile neutrinos are the leading candidate 
for warm dark matter and can be tested in x-ray searches from dark matter over densities, as the 
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center of the galaxy or nearby dwarf galaxies (talk by A. Kusenko). In the discussion session lead 
by P. Di Bari, various issues were discussed, in particular those related to the comparison between 
cosmological bounds and terrestrial ones, the estimation of systematic errors in the measurement of 
neutrino parameters from cosmological observations, and the impact of a non-standard evolution of 
the Early Universe. 



8.3. A theoretical perspective on lepton flavor physics at the TeV scale within i) 
SUSY models ii) extra-dimension models 

The origin of neutrino masses requires new particles and interactions which need to be embedded in a 
wider scheme to solve the hierarchy problem and explain the pattern of masses and mixing. Possible 
scenarios are super symmetric (talks by S. Lavignac, S. Khalil and A. Texeira), extra dimensions, and 
little Higgs models. One needs to consider how these models fit into a wider particle physics theory 
and if there are any model-independent signatures. These issues were discussed by A. Santamaria. 
Generically, the embedding of neutrino mass models in a wider scheme induces new signatures as LFV 
processes at low energy such as ^ — t- 67, fi — e conversion, and r — t- ej. 



8.4. Leptogenesis in the context of neutrino mass models: model dependent versus 
model independent considerations 

Leptogenesis (talk by A. Abada) is one of the favored mechanisms for the generation of the baryon 
asymmetry of the Universe and naturally takes place in see-saw models. In the presence of CP-violation 
the decays of heavy Majorana neutrinos in the Early Universe produce a lepton asymmetry which 
is partially washed-out and converted into a baryon asymmetry by sphaleron processes. Compelling 
questions arise (discussion session lead by E. Nardi): Under what conditions can a connection between 
low and high energy CP violation be established? Is it obvious to believe in leptogenesis if a low energy 
CP violation and lepton number violation is observed? In see-saw type I models, in general, there 
is not a direct connection between low energy parameters and high energy ones. This link arises if 
a theory of flavor is present, which explains the observed mixing pattern and reduces the number of 
free-parameters. For masses of the heaviest neutrino smaller than 10^^ GeV, flavor effects play an 
important role. In this case, generically a connection between leptogenesis and low energy CP- violating 
phases can be obtained in the sense that if CP-violation is found in future neutrino experiments a 
lepton asymmetry needs to be generated in the Early Universe. Leptogenesis can also take place in 
other see-saw models, for example see-saw type II models and in other scenarios of neutrino mass 
generation. 
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8.5. Interplay between neutrino masses and other phenomenological signatures 

Small neutrino masses require new physics beyond the Standard Model which can induce also other 
signatures at low energy. Additional sources of lepton flavor violation induce processes at low energy 
such as — )• 67, ^ — e conversion, r — )• 67 (talk by A. de Gouvea) at rates which could be close to 
present bounds. The LHC experiment will provide information on the physics at the TeV scale which 
might also be at the origin of neutrino masses (talk by T. Schwetz). Additional effects can arise in 
neutrino oscillations due to non-standard interactions of neutrinos with matter which can be tested in 
future long-baseline neutrino experiments (talk by E. Fernandez-Martinez). Neutrinoless double beta 
decay can be mediated not only by light Majorana neutrino masses but also by other mechanisms as 
heavy sterile neutrinos, supersymmetric R-parity violation. Once a positive signal is found, it will 
be critical to discriminate between the various mechanisms. S. Kom gave a contributed talk on LHC 
probes of SUSY neutrinoless double beta decay mechanism and C. Jackson on the issues related to 
understanding the mechanism of neutrinoless double beta decay with the SuperNEMO experiment. 
In the discussion session T. Ohlsson addressed the synergy and complementarity of these different 
experimental signatures in pinning down the mechanism at the origin of neutrino masses. 

8.6. Discussion on performance indicators in long baseline experiments 

Future long baseline neutrino oscillation experiments play a crucial role in determining the value of the 
unknown mixing angle ^13, the type of neutrino mass hierarchy and the existence of CP-violation in the 
leptonic sector. A rich experimental program is under consideration for the near future: conventional 
and superbeam experiments (MINOS, OPERA, T2K, NOvA, LBNE) are already taking data, under 
construction or at the R&D phase while even more ambitious projects such as neutrino factories and 
beta-beams are studied. The sensitivity of the various setups is studied in detail. T. Li discussed 
the physics reach of the low energy neutrino factory and J. Lopez Pavon the search for 3+1 sterile 
neutrinos in contributed talks. It is necessary to compare the reach of different experiments in an 
objective way and with comparable assumptions. These issues where discussed in the round table lead 
by K. Long. 

8.7. Conclusions 

With the discovery of neutrino oscillations, neutrino physics has opened a new window on the physics 
BSM. This workshop has given the opportunity to review the status of the studies in neutrino phe- 
nomenology, theory and astroparticle physics with particular focus on the theoretical motivations for 
future precise neutrino oscillation experiments. 
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